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Predicting and characterizing the crystal structure of materials is a key problem in materials re- 
search and development. We report the results of ab initio LDA/GGA computations for the following 
systems: AgAu, AgCd, AgMg, AgMo*, AgNa, AgNb*, AgPd, AgRh*, AgRu*, AgTc*, AgTi, AgY, 
AgZr, AlSc, AuCd, AuMo*, AuNb, AuPd, AuPt*, AuRh*, AuRu*, AuSc, AuTc*, AuTi, AuY, 
AuZr, CdMo*, CdNb*, CdPd, CdPt, CdRh, CdRu*, CdTc*, CdTi, CdY, CdZr, CrMg*, MoNb, 
MoPd, MoPt, MoRh, MoRu, MoTc*, MoTi, MoY*, MoZr, NbPd, NbPt, NbRh, NbRu, NbTc, 
NbY*, NbZr*, PdPt, PdRh*, PdRu*, PdTc, PdTi, PdY, PdZr, PtRh, PtRu, PtY, PtTc, PtTi, 
PtZr, RhRu, RhTc, RhTi, RhY, RhZr, RuTi, RuTc, RuY, RuZr, TcTi, TcY, TcZr, TiZr*, YZr* 
(* = systems in which the ab initio method predicts that no compounds are stable). A detailed 
comparison to experimental data confirms the high accuracy with which ab initio methods can pre- 
dict ground states. 
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I. INTRODUCTION 

First principles computation, whereby the properties 
of materials are predicted starting from the principles 
of quantum mechanics, is becoming well integrated with 
more traditional materials research. A list of ab initio 
studies on binary and ternary alloy phase stability up to 
1994 can be found in reference [1]. Since the earliest, 
completely ab initio computation of a binary phase di- 
agram [2], the approaches to compute the total energy 
of a solid have significantly improved, and computing re- 
sources have continued to become faster and less expen- 
sive. We believe that a point has been reached where, 
with a reasonable amount of resources, high throughput 
first principles studies of a large number of alloys can 
be performed [3-8]. In this paper we present the results 
of a first principles study of 14080 computed total ener- 
gies on 176 crystal structures in 80 binary alloys. To our 
knowledge this is the largest first principles study of its 
kind on alloys. As we have compared the results in every 
system to experimental compilations, this study also of- 
fers a statistical test on the accuracy of some current ab 
initio approaches in correctly predicting the structure of 
materials. 

For 89 compounds we find unambiguous agreement be- 
tween experiment and the ab initio computation (Table 
5), giving some indication of the predictive power of mod- 
ern ab initio electronic structure methods. For many sys- 
tems, verification of the ab initio results is difficult, as the 
systems have been poorly or incompletely characterized, 
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or only high-temperature information is available experi- 
mentally. For most of these system, we make predictions 
that are consistent with the limited available informa- 
tion. Even though our library of 176 crystal structures 
is, to our knowledge, the largest library of ab initio ener- 
gies ever produced, there are still 27 compounds for which 
we cannot verify the experimental structures as they are 
not in our library. We have not included such prototypes 
because they are extremely rare and complicated (many 
atoms per unit cell) (Tabic 8). 

Overall, we find remarkably few significant discrepan- 
cies between the ab initio predictions and the experimen- 
tal observations (Tabic 9). Based on the experimental 
data in references [9,10], we find only nine compounds 
for which LDA/GGA seems to predict the ground state 
incorrectly. For four of tlic^se nine systems, the experi- 
mental ground state is within less than lOmeV/atom of 
the ab initio ground state. For the remaining five sys- 
tems, there are at least two in which further investigation 
indicates that the experimental structure assignment is 
poorly justified, leaving three compounds in which a sig- 
nificant disagreement between experiment and ab initio 
LDA/GGA is likely. Such disagreements are addressed 
in Section (IV). 

We believe that the low ratio of unambiguous errors 
(3) to the number of unambiguous correct predictions 
(89) is encouraging, and establishes clearly the potential 
of predicting crystal structure correctly with ab initio 
methods. 

We also predict the stability of five new crystal struc- 
tures which, to the best of our knowledge, have not yet 
been observed in any system. An AB3 superstructure of 
the fee lattice, stable for CdPts, PdPta and PdaPt, an 
AB bcc superstructure for MoTi, an AB3 bee superstruc- 
ture for MoTig, MogTi, NbgTc, RuTia and TcTig, an 
A2B2 hep superstructure for RhRu, and an A2B4 hep su- 
perstructure for RhRu2 (Appendix (XI)). In addition, we 
find two new crystal structures which are not superstruc- 
tures of fee, bcc or hep: MoZrs and MosTi (for MosTi, 
MoZrs, and NbsRu) (Appendix (XII)). 



II. THE LIBRARY: ALLOYS AND STRUCTURES 

Binary alloys. Our calculated library contains 80 
binary intermetallic alloys. The alloys include the bi- 
naries that can be made from row 5 transition metals, 
as well as some systems with Aluminum, Gold, Mag- 
nesium, Platinum, Scandium, Sodium, Titanium, and 
Technetium. The alloys are: AgAu, AgCd, AgMg 
AgMo*, AgNa, AgNb*, AgPd, AgRh*, AgRu*, AgTc* 
AgTi, AgY, AgZr, AuCd, AuMo*, AuNb, AuPd, AuPt* 
AuRh*, AuRu*, AuSc, AuTc*, AuTi, AuY, AuZr, AlSc 
CdMo*, CdNb*, CdPd, CdPt, CdRh, CdRu*, CdTc* 
CdTi, CdY, CdZr, CrMg*, MoNb, MoPd, MoPt, MoRh 



MoRu, MoTc*, MoTi, MoY*, MoZr, NbPd, NbPt, 
NbRh, NbRu, NbTc, NbY*, NbZr*, PdPt, PdRh*, 
PdRu*, PdTc, PdTi, PdY, PtRh, PtRu, PtY, PtTc, 
PdZr, PtTi, PtZr, RhRu, RhTc, RhTi, RhY, RhZr, 
RuTi, RuTc, RuY, RuZr, TcTi, TcY, TcZr, TiZr*, YZr*, 
where the superscript * indicates the systems in which the 
high throughput ab initio method predicts that no com- 
pounds are stable: 57 alloys form compounds and 23 are 
non compound forming. 

Structures and their prototypes. The library con- 
tains 176 crystal structures. Many of these have the same 
structure type but with different compositions of occu- 
pancies, for example. AB3 and A3B (also AB and BA if 
the point groups of atomic positions of A and B are dif- 
ferent), so the number of distinct prototypes is 101. The 
various concentrations are listed in the table below. 
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TABLE 1. Compositions, concentrations and number of 
prototypes inside the library. The library has 176 struc- 
tures, and 101 distinct prototypes (* at composition AB, 
3 prototypes have different point groups in atomic posi- 
tions A and B, therefore they represent distinct structure 
types). 

Of such prototypes, 65 are chosen from the most com- 
mon intermetallic binary structures in the CRYSTMET 
database [11] and in the Pauling File [10]. Such proto- 
types can be described by their Strukturbericht designa- 
tion and/or natural prototypes [9,10]: Al, A2, A3, A4, 
A15, B,„ Bl, B2, B3, B4, B81, B82, BIO, Bll, B19, B27, 
B32, B33 (B/), Cc, C2, C6, CUb, C14, C15, C15b, C16, 
CIS, C22, C32, C33, C37, C38, C49, D0„, DO3, DO9, 
DOii, DO19, DO22, DO23, DO24, DI3, Dla, D2d, D73, D88. 
Llo, Lli, LI2, L60, Caln2, CuTe, CuZr2, GdSi2 (1.4), 
MoPt2, NbAs (NbP), NbPds, Ni2Si, O (with z=l/4), 
PU3AI (C03V), Ti3Cu4, W5Si3, YCd3, ZrSi2, 7-Ir. The 
rest of the structures (36) are fee, bcc or hep superstruc- 
tures. Twelve of these superstructures consist of stacking 
of pure A and B planes along some direction. Of such 
prototypes, 12 contain a stacking direction, therefore we 
can name them following the parent lattice and the stack- 
ing direction: 

-^[direction] 



LATTICE^ 



stacking 



(1) 



For example the designation FCC^^2B2 indicates a struc- 
ture prototype with FCC parent lattice and A2B2 
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stacking along direction [001]. The 12 prototypes of 
the hbrary that can be labeled in this manner are: 

^[011] pnnliii] 

FCC|^^3, FCCjI^jg. For the FCC superstructures, a con- 
version table between the lattice-stacking-direction, the 
Sanchez - de Fontaine notation [12], and the Lu et al. 
designations [13,14] is included below. 
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TABLE 2. Conversion table between the lattice- 
stacking-direction, the Sanchez - de Fontaine notation 
[12], and the Lu et al. designations [13,14]. 

The prototypes of stable structures, that cannot be 
described by Strukturbcricht designation, natural pro- 
totypes or lattice-stacking-direction convention, are de- 
scribed in Appendix XL The complete geometrical de- 
scription of all the fee, bee, and hep superstructures of 
the library can be found in reference [15]. 



III. HIGH-THROUGHPUT FIRST PRINCIPLES 
CALCULATIONS 

Ultra Soft Pseudopotentials LDA calculations 
(US-LDA). IMost of the energy calculations of the li- 
brary were performed using Density Functional Theory 
in the Local Density Approximation (LDA), with the 

Ceperley- Alder form for the correlation energy as param- 
eterized by Perdew-Zunger [18] with ultra soft Vander- 
bilt type pseudopotentials [19], as implemented in VASP 
[20]. Calculations arc at zero temperature and pres- 
sure, and without zero-point motion. The energy cut- 
off in an alloy was set to 1.5 times the larger of the 
suggested energy cutoffs of the pseudopotentials of the 
elements of the alloy (suggested energy cutoffs are de- 
rived by the method described in [20]). Brillouin zone 
integrations were done using 2000/ (number of atoms in 
unit cell) k-points distributed as uniformly as possible 
on a Monkhorst-Pack mesh [21,22]. For a sample set of 
calculations, we verified that with these energy cutoffs 
and k-points meshes the absolute energy is converged 



to bettor than 10 meV/atom. Energy differences be- 
tween structures are expected to be converged to much 
smaller tolerances. Spin polarization was used only for 
magnetic alloys. For the non-magnetic ones, spin polar- 
ization was used only to verify the energies of the ground 
state structures. All structures were fully relaxed. With 
such methodology, we consider degenerate structures the 
ones that have energies closer than about 5meV/atom. 

We have decided to calculate the entire library within 
the LDA formalism instead of the PAW-GGA method 
(discussed below) because of the faster execution speed 
of LDA. 

PAW-GGA calculations. When several structures 
are in close competition for the ground state, we also per- 
formed calculations in the Generalized Gradient Approxi- 
mation (GGA), with Projector Augmented- Wave (PAW) 
pseudopotentials, as implemented in VASP [20,23,24]. 
One would expect the PAW-GGA approach to be some- 
what more accurate than the US-LDA. For the GGA cor- 
relation energy, we used the Perdew-Wang parameteriza- 
tion (GGA-PW91) [25]. Similar to the US-LDA case, all 
PAW-GGA calculations are performed at zero tempera- 
ture and pressure, and without zero-point motion. The 
energy cutoff in each calculation was set to 1.75 times the 
larger of the suggested energy cutoffs of the pseudopoten- 
tials of the elements of the alloy. Brillouin zone integra- 
tions were done using at least ~3500/ (number of atoms 
in unit cell) k-points distributed as uniformly as possible 
on a Monkhorst-Pack mesh [21,22]. Spin polarization was 
used in all calculations. We expect such calculations to be 
able to distinguish energies of structures previously de- 
generate. This is because the increased density of the k- 
points mesh, the reduced radial cutoff of the PAW poten- 
tials versus the ultrasoft pseudopotentials, and the GGA 
correlation energy implementation. With such method- 
ology, we consider degenerate structures the ones that 
have energies closer than <lmcV/atom. 

For some structures both the US-LDA and PAW-GGA 
calculations are presented. To avoid confusion, we specify 
the type of calculation in brackets: "(us-lda)" or "(paw- 
gga)". In addition, all the results oj the PAW-GGA cal- 
culations are described in italics. 

Symmetries of the pure elements. Both the US- 
LDA and PAW-GGA calculations reproduce the correct 
experimental crystal structures of the pure elements at 
room temperature. The only exception is Sodium. In the 
two formalisms, Na-hcp is very slightly favored over Na- 
bcc and Na-fcc, in agreement with other first principle 
calculations [26 30]. At room temperature Sodium has 
the bcc structure, and undergoes a martensitic transfor- 
mation below 35K, to a closed packed structure [31-38]. 
Therefore our results, Na-hcp stable and the very small 
energy differences with Na-bcc and Na-fcc, are consistent 
with the behavior at low temperature only. 
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High-throughput computing. To perform the high 
number of calculations, we have implemented a set of au- 
tomatic tools to prepare initial data, perform calculations 
and analyze results. 

Preparation of input files. All input files are prepared 
by starting from the templates based on the prototype 
chemistry. The volume of the unit cells is determined as 
a linear combination of atomic volumes of pure elements 
with Vegard's law [39] . The internal position of the atoms 
have been taken from geometrical configurations (for fee, 
bcc, hep superstructures) or have been extracted from 
the CRYSTMET and the Pauling File databases [11,10]. 

Ab initio Calculations. The calculations are performed 
by a high-throughput tool, called AFLOW (automatic 
flow), which searches inside the library for potential in- 
put files and performs the proper ab initio calculations. 
AFLOW is able to balance the CPU loads in a multi- 
processor and cluster environment, maximizing the to- 
tal throughput of the process [15]. AFLOW extracts the 
energy and initiates further relaxations with the atomic 
positions and unit cell geometries obtained from VASP. 
Therefore, all the structures are relaxed at least twice, 
and many (~ 10% of the total) are relaxed three or more 
times, until proper geometrical and energetic convergence 
is obtained. To conclude, for a library of 176 crystal 
structures in 80 binary alloys, with a total of 14080 struc- 
tures, AFLOW performed 32,402 VASP calculations. The 
total computing time used for this project was ^ 9.05-10^ 
CPU seconds (~ 28.7 CPU years), which was spread over 
a large set of computers [40] . Given the different types 
of CPUs used, we can only estimate the total amount of 
computation for the project at roughly ~ 1.2 • 10^ TER- 
AFLOP. 

Collection and error check of the results. Once each 
structure has been calculated, space group symmetry, 
bond distances, coordination numbers and unit volumes 
are compared for all the structures of that alloy with 
the same concentration in the library. Frequently, sev- 
eral structures, starting from different prototypes, relax 
to the same final structure with the same energy. There- 
fore, to identify the correct stable phase, it is manda- 
tory to be able to determine the final relaxed structure 
with the highest possible reliability. Details of the high- 
throughput error checking tool can be found in reference 
[15]. A few stable structures remain unidentified. For 
those we report the unit cell and atomic positions in Ap- 
pendix XI. Such prototypes might be new phases to be 
checked experimentally. 

Calculation of the formation energies and the convex 
hull. The formation energy for each structure is deter- 
mined with respect to the most stable structure of the 
pure elements. To determine the ground states of a sys- 
tem one needs to find, as a function of composition, the 
ordered compounds that have an energy lower than any 



other structure or any linear combination of structures 
that gives the proper composition. This set of ground 
state structures forms a convex hull, as all other struc- 
tures have an energy that falls above the set of tie lines 
that connects the energy of the ground states. In ther- 
modynamical terms, the convex hull represents the Gibbs 
free energy of the alloy at zero temperature. 

Comparison with experiments. To compare the 
ab initio results to experimental information, we relied 
on the information and references in the Binary Alloy 
Phase Diagrams (Massalski [9]) and in the Pauling File 
[10]. Though in some cases references not included in 
these were also used, no systematic approach to go be- 
yond information in the Pauling File or Binary Alloy 
Phase Diagrams was used. 

IV. DISCUSSION AND SUMMARY OF RESULTS 

In comparing the stable structures predicted by the 
ab initio computations with available experimental in- 
formation, we have attempted to classify the results in 
a few distinct categories. Table 5 gives the compounds 
where the ab initio result and experiments are in un- 
ambiguous agreement. The fact that there are a large 
number of compounds (89) in Table 5 is a positive state- 
ment about the accuracy of LDA/GGA in capturing the 
close energetic competition between the 176 structures in 
our library. For one of the systems in Table 5 (PdTia) 
we took the liberty of modifying the experimental re- 
sult [9] which shows the A15 structure (stoichiometry 
AB3) as a line compound at composition PdTi4. While 
off-stoichiometric compounds are obviously possible, this 
usually goes together with significant width of the single- 
phase field. Hence, we concluded (maybe erroneously) 
that the placement of A15 at composition PdTi4 in ref- 
erence [9] is likely a typographical error. 

Table 6 shows compositions which are experimentally 
known to form compounds, but lack a complete identi- 
fication of the structure type. Hence, the ab initio pre- 
diction should be seen as a likely crystal structure for 
the compound. In most cases, the ab initio structure is 
consistent with the conditions imposed by the limited ex- 
perimental data (e.g. bcc-ordering is seen for Ag-Cd, and 
we predict B19; or C37 is speculated for Au-Sc, in agree- 
ment with what we predict). Hence, it is likely that with 
further experimental characterization, many of these sys- 
tems would move to Table 5. 

Table 7a-b contains compositions which are experimen- 
tally characterized as solid solutions, high-temperature 
two-phase regions, or have not been studied at all in a 
particular composition range. Because of the lack of low- 
temperature information in these systems, an unambigu- 
ous statement about the accuracy of the LDA/GGA pre- 
diction cannot be made. For some systems, the ab initio 
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result is very consistent with the experimentally observed 
behavior. For example, in Ag-Au we find many fee or- 
dered superstructures with low formation energies, indi- 
cating a weak ordering interaction. So it is likely that the 
ab initio predicted structures are stable at low tempera- 
ture, but disorder into a solid solution at elevated tem- 
perature, a rather common occurrence for noble-metal 
alloys. For other systems in this table the comparison be- 
tween experiments and ab initio is more troubling. For 
example, while Mo-Ti is experimentally described as a 
two-phase system with no known compounds down to 
400°C, we find very strong compound formation. While 
it is possible that these compounds disappear through 
peritectoid reactions below 400° C, the large negative for- 
mation energies make this unlikely. It is more likely that 
either the experimental characterization of this system, 
or the ab initio result is significantly in error. A system 
such as Mo-Nb, on the other hand, is not at all studied 
experimentally, and no statement about the accuracy of 
our prediction in this system can be made. 

Table 8 contains the compounds whose crystal struc- 
ture is not present in the library, and hence, it is not 
possible to make a comparison with the ab initio data. 
In some cases (e.g. AuioZry), this is because the crystal 
structure is uncommon and has a large unit cell, making 
it less worthwhile to include it in the library. In other 
cases (e.g. RhsTi) the structure is simply unknown. Ta- 
ble 8 also includes structures that are only stable at off- 
stoichiometric compositions (e.g. B19 in Au-Ti), as we 
did not include disorder in the library prototypes. 

Finally, Table 9 contains the compositions for which 
there seems to be a clear disagreement between the ex- 
perimental data in references [9,10] and the ab initio re- 
sults. Because these systems can point at cither experi- 
mental errors, or shortcomings of the ab initio method, 
we discuss them here in more detail. Note that ab initio 
errors are most likely due to the LDA or GGA approxi- 
mations or the pseudopotentials being used. In general, 
one should consider that our calculations only produce 
the zero-K energy, whereas all the experimental results 
are at non-zero temperature. Given that entropic differ- 
ences between structures can be of the order 0.1 ~ l.Ofcs 
per atom [381], very small energy differences between 
the experimentally observed structure and our ab ini- 
tio results, could be reversed at elevated temperature. 
In particular, AuNbs-AlS, NbRh3-Ll2, NbRu3-Ll2 and 
NbRu'-Llo may be in this category. This phenomenon 
is quite common. As an example, Wolverton and Ozilins 
have shown that, in the case of the aluminum-copper 
system, the vibrational entropy difference is responsible 
for stabilizing the 6 phase AI2CU (tetragonal C16) over 
the competing AI2CU-6'' phase (distortion of 6'c-Cl, the 
cubic fluorite phase with CaF2 prototype), which has 
the lowest energy, and, therefore, is the stable struc- 



ture at zero Kelvin [382]. At temperatures higher than 
T~ 150 - 200°C, AI2CU-6' is more favorable. 

A. Discussion of disagreement between experimental and 
ab initio results 

AuNbs. The structure AuNbs-AlS has been re- 
ported a munbcr of times [9,10,127 141]. In particular, 
Roschcl et al. [132] observed the A15 structure to grow 
in quenched samples. Ball-milling samples of the mate- 
rial to form disordered bee, revert upon heat treatment 
in an exothermic reaction to A15, suggesting that A15 
could be the ground state for this composition at low 
temperature;. The relatively small energy difference be- 
tween A15 and the AuNb2^Nb tie-line in our calculation 
could lead to finite temperature stability of A15. Alter- 
natively, the ^7meV/atom energy difference (US-LDA 
and PAW-GGA) could simply be an error of the ab ini- 
tio approach. 

AuY. While we predict the B33 structure to be sta- 
ble by 26meV/atom below the B2 structure, this system 
is experimentally listed as having a B2 structure. This 
experimental classification is based on a paper by Chao 
et al. [170,10]. Chao's work mentioned that they only 
see B2 when performing a rapid quench. Further work 
by Kusma [172] and Dwight [158] did not observe the B2 
structure. The experimental results in this system are 
therefore somewhat suspect. It is possible that B2 is a 
high temperature structure while B33 is the low temper- 
ature form. However, it should also be pointed out that 
the mixing energies in this system are very large ( ~leV), 
and that the potential error of ~25meV/atom (US-LDA 
and PAW-GGA) is a small fraction of this. 

CdsNb. The structure of this compound is one of 
the more interesting discrepancies between ab initio and 
experiments to emerge from our study. Based on ex- 
perimental work of Von HoUeck [178], a single LI2 com- 
pound is claimed to be stable in this system. Com- 
putationally, no compounds are stable at all, and we 
find this system to be immiscible. The compound with 
lowest energy is a bcc superstructure at composition 
CdNb, and energy '~55meV/atom above the Cd^Nb tie- 
line ('^62meV/atom with PAW-GGA). The CdsNb-Lla is 
70meV/atom above that same tie-line (> lOOmeV/atom 
with PAW-GGA). It is not common for LDA/GGA to 
make such large and qualitative errors in metals, and fur- 
ther investigation of this system is suggested, as it will 
likely lead to a reconsideration of the experimental clas- 
sification, or to some novel understanding of the errors 
that ah initio methods can produce in metals. 

CdPts. Experimentally this compound is claimed to 
form in the LI2 structure, while we predict it to form 
in a new structure for which no prototype yet exists 
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(CdPtf °*° in Appendix VIII). We believe the experimen- 
tal assignment to be uncertain. It is based on a single 
paper by Nowotny [181] which actually does not observe 
the LI2 structure in this system, but assigns it by chem- 
ical similarity with the PtaZn system. 

NbRha. This compound is characterized as an 
LI2 structure, whereas our ab initio results find the 
AI3PU prototype to be 8meV/atom lower in energy 
(5.3meV/atom with PAW-GGA). This small energy dif- 
ference could very well be an LDA/GGA error, or could 
be reversed at elevated temperature by entropic effects. 

NbRua. NbRus is similar to the discrepancy of com- 
pound NbRhs. Experimentally it is foimd to be an 
LI2 structure, whereas we find a DO24 structure that is 
SmeV/atom below the LI2. With PAW pseudopotentials 
and the GGA exchange correlation corrections the DO24- 
LI2 energy difference reduces to 2.5meV/atom. Because 
this number is so small, extremely large k-points sets and 
high energy cutoff were used to converge it. The results 
indicate that LI2 and DO24 are all but degenerate in this 
system. Since DO24 and LI2 are very similar (DO24 is a 
periodically anti-phased version of LI2) subtle entropic 
effects can easily reverse the stability between these two. 

NbRu. Experimentally this system is classified to 
have the LIq structure, whereas the ab initio results 
indicate an unusual two-phase region between NbsRu 
(DO3) and NbRu2 (C37) covering this composition 
[9,10,233,237]. The most recent experimental result is 
based on the assignment by Chen [237] who observes a 
disordered bcc solution at high temperature and a trans- 
formation to a tetragonal structure near 1000°C. Be- 
tween 720~920°G they see a further symmetry reduction 
to an orthorhombic phase. Hence, even if the tetrago- 
nal structure below 1000°G would be Llo (which is only 
speculated by Chen) its further transformation to an or- 
thorhombic structure indicates that LIq is not the ground 
state structure. Based on a diffusion couple study. Hurley 
[235] believes that a two-phase region exists near 50% Ru, 
though the results seem to be irreproducible. Besides the 
difficulty with interpreting the experimental data, signif- 
icant problems also exists on the ab initio side with this 
compound. The PAW-GGA result is substantially dif- 
ferent from the US-LDA result. GGA gives Llo only 
4meV/atom above the tie-line (versus 20meV/atom in 
LDA), whereas B19 seems to be >100meV/atom above 
the tie-line (versus 13mcV/atom with LDA). 

Pt Y. Experimental assignment of the structure of this 
compound as B27 is based on two papers [277,316]. The 
ab initio work finds B33 to be lowest in energy with B27 
higher by 50meV/atom (60meV/atom with PAW-GGA). 
Given that the experimental results seem sound, this is 
probably a true ab initio error. 

PtaZrs. The work of Schubert et al. [323] establishes 
Pt3Zr5 as a DBg structure, whereas our ab initio (LDA) 



results put the DSs structure 36mcV/atom above the tie- 
line formed by PtZr2-C16 and PtZr-B33. The GGA re- 
sult is 26meV/atom above the tie-line. In each case the 
WsSia structure is lower in energy than the DSg. Given 
the large mixing energies in this system (of the order of 
1 eV/atom), these ab initio errors seem plausible. 

V. CONCLUSIONS 

Overall, the comparison between experimental data 
and ab initio results is encouraging. A large number of 
the ground states are predicted correctly, even though 
significant competition exists between various structures, 
indicating that relative energy differences are well repro- 
duced by LDA/GGA. In many cases, the direct compari- 
son between ab initio and experiments is difficult, as sys- 
tems have often not been studied completely, or have not 
been studied to low enough temperature to make a re- 
liable statement about the ground state structures. For 
only a small number of cases, there is a clear discrep- 
ancy between ab initio and experiments (Table 9). On 
further investigations, we find that some of the experi- 
mental assignments are poorly justified (e.g. CdPts). In 
a few of the "error systems", the energy difference be- 
tween the observed structure and the ab initio structure 
is small, which makes it fall within the uncertainty one 
would expect from LDA/GGA or from a lack of entropic 
considerations. 

A few systems display, in our opinion, such a signif- 
icant discrepancy between available experimental data 
and ab initio results as to make them worthwhile for fur- 
ther study. These include CdsNb (LI2 observed, but no 
compounds from ah initio), Nb-Ru and PtaZrs (respec- 
tively Llo and DSs observed, but no stable compounds 
at that composition in ab initio), PtY (significant error 
between B33 and B27) and the complete Mo-Ti system 
(which from experiments appears to be a miscibility gap, 
but in the ab initio results display a number of novel and 
strongly stabilized compounds). These systems should 
be investigated further, cither to find experimental er- 
rors, or to shed new light on ab initio problems. 

Our findings are summarized in Table 3. By compar- 
ing all the available ab initio calculations with the ex- 
perimental results [9,10], we obtain 284 potential struc- 
ture comparisons. A subset of these, (48), are not avail- 
able due by the lack of the proper prototype or the un- 
known experimental compound. The rest of the struc- 
ture comparisons, (236), can be divided in 89 agreements 
in ordered systems (Nac), 21 agreements in immiscible 
systems (Nai), 9 disagreements (Nj), 21 good predictions 
(Ngp), and 96 possible predictions (Npp). This division is 
summarized in the following table. 
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Total Experiments comparisons: 284 


Available comparisons: 236 




Experimental 


Ab initio 




# 


Table 


Agreements 


Compound 


Compound 


Nac 


89 


5 


Agreements 


Immiscibility 


Immiscibility 




21 


3 


Disagreements 


Comp./ Immisc. 


Incorrect 


Na 


9 


9 


Good 
predictions 


Comp. unknown, 
uncertain, guessed 


Compound 
prediction 


Ngp 


21 


6 


Possible 


not-studied. 


Possible 


Npp 


96 


7a-b 


predictions 


solid-solution, 
two-phase region 


compound 






Unavailable comparisons: 48 


Reason 


# 


Table 


unavailable ab 


initio prototype 






21 


8 


unknown experimental compound 


27 


8 



TABLE 3. Summary of the 284 structure comparisons 
between experimental and ab initio calculations. When 
available, we use the PAW-GGA calculations, otherwise 
we use the US-LDA ones. If only US-LDA calculations 
were used, we would obtain N^J-'^=84: agreements and 
iVj"^^=14 disagreements. 



The reliability of our ab initio method can be defined 
as the fraction of correct compounds (agreements) over 
the number of accessible compounds (agreements and dis- 
agreements): 

Such quantity, rjc, contains the reliability of the experi- 
mental results, which can be removed by considering only 
the true unambiguous disagreements oi Table 9 {N^ = 3). 
Hence the reliability of the method becomes: 

= « 96.7%. (3) 

Such quantities, rjc and rj*, can be considered as the prob- 
abilities that the most stable ab initio compounds repro- 
duce the correct experimental compound. The quantity 
T]* should be used if the experimental compound is cer- 
tain. 

If we extend the argument to the agreements with ex- 
perimental compounds and immiscible systems, then the 
reliabilities become: 

where the quantity r]* should be used if the experimental 
compound or immiscibility is certain. Such quantities, 
T] and 77*, can be considered as the probabilities that ab 
initio results reproduce the correct experimental com- 
pounds or immiscibilitics. 

Given the large number of calculations of this project, 
we believe r], 77*, 77c, and 77* to be good estimations of 
state of the art pscudopotential ab initio accuracies in 
phase stability prediction. 
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VI. ALLOYS WITHOUT AB INITIO COMPOUNDS 

Table 4 gives the alloys for which we find no com- 
pounds with negative formation energy, and the structure 
with lowest formation energy in the system. All of these 
agree with experiments except for the ones described be- 
low. 



System 


Structure 
with lowest Ey 


Ef (us-lda) 
mcV /atom 


References 
[9-11] and 


Ag-Mo 


AgMo3-BCC 
proto. #72 in [16] 


208 


[81,82] 


Ag-Nb 


AgNbs-BCC 
proto. #72 in [16] 


62 


[86,87] 


Ag-Rh 


Till! 

AgRh3-FCC[i'„'l 


116 


[50,94-96] 


Ag-Ru 


AgRu3-HEX 
proto. #127 in [16] 


199 


[41,43,96-98] 


Ag-Tc 


AgTc-HEX 
proto. #126 in [16] 


147 


[99] 


Au-Mo 


AuMo3-D0i9 


75 


[41,46,118-122] 


Au-Pt 


Au2Pt-FCC'{Rj 


9 


[41,117,119,146-148] 


Au-Rh 


AuRh3-FCC'°R] 


91 


[116,119,149,150] 


Au-Ru 


AuRu3-HEX 
proto. #125 in [16] 


165 


[116,119,151,152] 


Au-Tc 


AuTc-HEX 
proto. #126 in [16] 


74 


[116,159] 


Cd-Mo 


CdsMo-HEX 
proto. #128 in [16] 


173 


[9] 


Cd-Nb 


CdNb-BCC 
proto. #71 in [16] 


55 (us-lda) 
58 (paw-gga) 


[178] 

see note "Cd-Nb" 


Cd-Ru 


CdRu-FCC^j;^'^, 


88 


see note "Cd-Ru" 


Cd-Tc 


CdTc-Bll 


92 


[182] 


Cr-Mg 


CrMg-Bll 


251 


[45,191] 


Mo-Tc 


M05TC-AI5 


53 


[44,217] 
see note "Mo-Tc" 


Mo-Y 


M02Y-CI5 


40 


[222,223] 


Nb-Y 


Nb2Y-C49 


132 


[244-248] 


Nb-Zr 


NbZr2-HEX 
proto. #129 in [16] 


22 


[249-253,90] 
[238-240,254-258] 


Pd-Rh 


PdRh3-FCC[{^':| 


35 


[13,263-270] 


Pd-Ru 


PdRu-HEX 

proto. #126 in [16] 


44 


[271-274] 


Ti-Zr 


Ti2Zr-B82 


18 


[48,49,371-374 


Y-Zr 


YZrs-HEX 
proto. #141 in [16] 


35 


[49,375-380 



TABLE 4. Systems without intermetallic compounds. 
The second and third columns give, for each system, the 
structure with the lowest formation energy E f (US-LDA 
calculations). The table contains 23 entries. 21 of these 
are in agreement with experiments. References include 
both experimental and ab initio studies. 

Cd-Nb. One compound, Cd3Nb-Ll2, has been re- 
ported for the system Cd-Nb [9,178]. However, we did 
not find any stable phase (Cd3Nb-Ll2 has formation en- 
ergy of ~70meV/atom). The disagreement of compound 
Cd3Nb-Ll2 is further discussed in the Section (IV). 

Cd-Ru. The authors are not aware of any experimen- 
tal result for the Cd-Ru system. 



Mo-Tc. Two compounds have been suggested 
based on thermodynamic models: /3]VIo2Tc3-A15 and 
(j]VIo3Tc7-D8f, [9,44]. However, we did not find any sta- 
ble phase for the system, nor can we check the suggested 
compounds, since our library does not have the a phase 
or the off-stoichiometry A15. 
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VII. ALLOYS WITH AB INITIO COMPOUNDS 

Ag-Au (Silver - Gold). The system Ag-Au has 
not been studied in great detail at low temperatures, 
and no intermetallic compounds have been reported 
[9,10,52-55,119]. The solid is reported to be short- 
range ordered fee, though it is suggested that long- 
order might exist at low temperature. At low temper- 
ature we find several stable compounds: Ag4Au, AgaAu, 
Ag2Au, AgAu-Llo, AgAu2 and AgAua. In our electronic 
structure approach, the ground states are degenerate for 
AgsAu, Ag2Au, AgAu2 and AgAus. Our best candidates 
are LI2, DO23, AlgPu, NbPdg, DO22 and DO24 for AggAu, 
C37 and MoPt2 for Ag2Au and AgAu2, and LI2, DO22 
and DO23 for AgAus, as shown in figure (1). The consid- 
erable degeneracy of the ground states indicates that only 
very small effective interactions exist between Ag and Ag, 
consistent with the experimentally observed complete 
solid solubility. For composition Ag4Au, the structure 
Dla is degenerate with the tie-line of the two-phase re- 
gion Ag<->Ag3Au. We conclude that further experimental 
and theoretical investigations are necessary to determine 
the behavior of AgAu. To address the degenerate struc- 
tures, we further investigate Ag^Au, Ag2Au, AgAu2 and 
AgAu^ with PAW-GGA potentials, as described in Sec- 
tion (III). With PAW, for composition Ag^Au, LI2 is the 
most stable compound and DO22, DO23, Al^Pu, NbPd^, 
DO24 are higher by 0.9, 1.6, 2.2, 2.4, and S.OmeV/atom, 
respectively. For both composition Ag2Au and AgAu2, 
C37 is the most stable compound and MoPt2 is higher 
by 1.9meV/atom and 4-3meV/atom, respectively. For 
composition AgAu^, LI2 is the most stable compound 
and DO23 md DO22 are higher by 3.4meV/atom and 
3.9meV/atom, respectively. 

Other ab initio studies, relevant for this system, can 
be found in references [56-66]. 



Ag-Au system 


Low Temperature Phases comparison ehart 


% Au 


(Massalski [9]) 


Ah "JTi '?/') ^^ 

result 


20 


short-range order 
> 950° C 


Ag4Au-Dla/tie-line 




short-r&iigG order 
> 950° C 


A O-n A 11 

LI2/DO23/AI3PU/ (us-lda) 

LI2 stable (paw-gga) 
D022^0.9meV/at. 
D023^1.6meV/at. 

Al3Pvr^2.2meV/at. 

NbPd3r^2.4meV/at. 
D024^3.0meV/at. 

above LI2 (paw-gga). 


33.3 


short-range order 
> 950° C 


Ag2Au-C37/MoPt2 (us-lda) 
C57 stable (paw-gga) 
Mo Pt2 ~ 1 . 9me V/at. 
above C37 (paw-gga). 


50 


short-range order 
> 950° C 


AgAu-Llo 


66.6 


short-range order 
> 950° C 


AgAu2-C37/MoPt2 (us-lda) 
C37 stable (paw-gga) 
MoPt2^4.3meV/at. 
above C37 (paw-gga). 


75.0 


short-range order 
> 950° C 


AgAu3 
LI2/DO22/DO23 (us-lda) 
LI2 stable (paw-gga) 
D023^3.4meV/at. 
D022 ^3.9meV/at. 
above LI2 (paw-gga). 



0.01 




FCC 



40 60 
Atomic Percent Gold 



100 
Au 



FIG. 1. AgAu (Silver - Gold) ground state convex hull. 
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Ag-Cd (Silver - Cadmium). At composition 
AgCds we find DOig stable, while only disordered hep 
solid solution has been reported at this composition. 
Hence, DOig, an ordered hep superstructure, probably 
is stable at low temperature in this system. At compo- 
sition AgCd2 we find a previously unknown compound 
AgCd2-ZrSi2. The low temperature near AgCd sto- 
ichiometry, is reported to be ordered bcc in Massalski 
[9]. In particular, the Pauling File reports AgCd-B2 
[10,67-70]. We find B2, B19 and B27 to have degenerate 
energies (B19 is quite common for Cd alloys). The ex- 
perimental phase diagram only displays solid solution on 
the Ag-rich side [9,10]. We find C37 at Ag2Cd and a sta- 
ble phase with stoichiometry AgaCd. The ground state 
of AgsCd is degenerate: three distinct structure (DOig, 
DO22, DO24) have similar energy. To address the degen- 
erate structures, we further investigate AgCd and Ag^ Cd 
with PAW-GGA potentials, as described in Section (HI). 
With PAW, for compound AgCd, B19 is the most stable 
structure, and B27 and B2 are higher by 2.8meV/atom 
and S.J^meV/atom, respectively. For compound Ag^Cd, 
DO22 and DO24 are still degenerate ground states and 
DOig is higher by only 2.2meV/atom. It is possible that 
these compounds have not yet been observed, or that they 
have low ordering temperatures. 



Ag-Cd system 


Low Temperature Phases comparison chart 


Composition 
% Cd 


Experimental 
(Massalski [9]) 


Ah initio 
result 


25.0 


fee solid solution 


AggCd 
DO19/DO22/DO24 (us- Ida) 
DO22/DO24 (paw-gga) 
D0i9r^2.2meV/at. 
above DO22/DO24 (paw-gga) 


33.3 


fee solid solution 


Ag2Cd-C37 


48.5 to 50 


P'-hcc 
B2 [10] 


AgCd-B2/B19/B27 (us-lda). 
B19 stable (paw-gga) 
B27^2.8meV/at. 
B2 r^3.4meV/at. 
above B19 (paw-gga). 


66.6 


none 


AgCd2-ZrSi2. 


64.5 to 81 


hep solid solution 


AgCd3-D0i9 



-0.1 




Ag-Mg (Silver - Magnesium). The phase diagram 
of the system Ag-Mg is known with reasonable accu- 
racy [9,10,71,72,75,77]. Our ab initio method confirms 
the stability of Ag]V[g-B2. In the Ag-rich side of the 
phase diagram, at stoichiometry AgaMg, we find DO23 
and DO24 to be degenerate. While Massalski [9] indi- 
cates AgaMg to be LI2, it is known from detailed High 
Resolution Electronic Microscopy (HREM) that AgaMg 
forms Long Period Superstructure (LPS) modulations of 
LI2 [73,74] (LPS DO23 is present in our library). For a 
detailed discussion of this system see Kulik et-al. [76]. 
LI2 is 15.2meV/atom above DO23/DO24. On the Mg-rich 
side, our library does not have prototypes at 80% Mg 
concentration, so we are not able to find any hexagonal 
AgMg4 (with unknown prototype [75,77]). However we 
find a stable compound AgMga. Two phases, AgMga- 
DOa and AgMga-DOig, have degenerate energies. This 
is important, since experiments suggest that AgMg4 was 
probably identified as hexagonal AgMga in early stud- 
ies. Unfortunately, we cannot determine the correct ab 
initio solution since we lack the hexagonal AgMg4 proto- 
types in our library. To address the degenerate structures, 
we further investigate Ag^Mg and AgMg^ with PAW- 
GGA potentials, as described in Section (III). For com- 
pound Ag^Mg, with PAW, DO23 is the most stable com- 
pound (in agreement with [73,74,76]), -DO24 and LI2 are 
l.lmeV/atom and ll.SmeV/atom higher than DO23, re- 
spectively. In addition, for compound AgMg^, DOig is the 
most stable structure and DOa is '^6.8meV/atom higher 
than DOig. 

Other ab initio studies, relevant for this system, can 
be found in references [78-80]. 



Ag-Mg system 


Low Temperature Phases comparison chart 


Composition 
% Mg 


Experimental 
(Massalski [9]) 


Ab initio 
result 


25.0 


LI2 

DO23 [73,74,76] 


Ag3Mg-D023/D024 (us-lda) 
Ll2~15meV/atom above. 
£•023 stable (paw-gga) 
DO2A ^l.lmeV/at. 
LI2 r~.ll.3meV/at. 
above DO23 (po-w-gga) 


35.5 to 65.4 


B2 


AgMg-B2 


75.8 to 78.2 


cF* (unknown) 


AgMg3 
DO19 (hP8)/D0a (0P8) (us) 
DOig stable (paw-gga) 

DOa ^6.8meV/at. 
above DOig (po-w-gga) 


80.0 


hP* (unknown) 


unavailable (see text) 



Ag 



20 40 60 80 

Atomic Percent Cadmium 



100 
Cd 



FIG. 2. AgCd (Silver - Cadmium) ground state convex hull. 
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B2 



20 40 60 80 100 

Ag Atomic Percent Magnesium IVlg 

FIG. 3. AgMg (Silver - Magnesium) ground state convex hull. 



Ag-Na (Silver - Sodium). The phase diagram of the 
system Ag-Na is known with reasonable accuracy and has 
only one intermetallic compound [9,10,83-85]. We con- 
firm the stability of Ag2Na-C15 and find almost all the 
other compounds to have positive formation energies en- 
ergies. Both with US-LDA and PAW-GGA, Na-hcp is 
very slighly favored over Na-bcc and Na-fcc, in agree- 
ment with other first principle calculations [26-30]. At 
room temperature Sodium has the bcc structure, and 
undergoes a martensitic transformation below 35K, to a 
closed packed structure [31-38]. Therefore our results, 
Na-hcp stable and the very small energy differences with 
Na-bcc and Na-fcc, are consistent with the behaviour at 
low temperature. 



Ag-Na system 


Low Temperature Phases comparison chart 


Composition 
% Na 


Experimental 
(Massalski [9]) 


Ab initio 
result 


33.0 


C15 


Ag2Na-C15 




20 40 60 80 100 

Ag Atomic Percent Sodium Na 

FIG. 4. AgNa (Silver - Sodium) ground state convex hull. 
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Ag-Pd (Silver - Palladium). The system Ag-Pd 
has not been studied in great detail and no intermetalhc 
compounds have been reported [9,10,41,88]. The sohd is 
reported to be disordered fee. At low temperature we find 
three stable compounds: AgPd-Lli, Ag2Pd and AgsPd. 
In our computations, for Ag2Pd and AgsPd, the ground 
states are degenerate: our best candidates are C49 or C37 
for Ag2Pd, and LI2 or DO22 for AgaPd, as shown in figure 
(5). To address the degenerate structures, we further in- 
vestigate Ag2Pd and Ag^Pd with PAW-GGA potentials, 
as described in Section (III). With PAW, for composi- 
tion Ag2Pd, C37 is the most stable compound and C49 
is higher by ^meV/atom. For composition Ag^Pd, LI2 
and DO22 remain degenerate. In fact, DO22 has an en- 
ergy O.^meV/atom greater than LI2, too small compared 
to the numerical accuracy of the ab initio calculation. 

Other ab initio studies, relevant for this system, can 
be found in references [89-93]. 



Ag-Pd system 


Low Temperature Phases comparison chart 


Composition 
% Ag 


Experimental 
(Massalski [9]) 


Ah initio 
result 


50 


solid solution 
> 900° C 


AgPd-Lli 


66.6 


solid solution 
> 900° C 


AgzPd 

C37/C49 (us-lda) 
C37 stable (paw-gga) 

C49r^4meV/at. 
above C37 (paw-gga) 


75.0 


solid solution 
> 900° C 


AgsPd 
LI2/DO22 (us & paw-gga) 



Ag-Ti (Silver - Titanium). The stability of AgTi2- 
Cllfe and AgTi-Bll is confirmed [9,10,100,101], and no 
other stable phases are found computationally. Hence, 
we conclude that the low temperature part of the phase 
diagram of AgTi is probably accurate. 



Ag-Ti system 


Low Temperature Phases comparison chart 


Composition 

% Ag 


Experimental 
(Massalski [9]) 


Ab initio 
result 


33.3 




AgTi2-Cll6 


50 


Bll 


AgTi-Bll (7CuTi) 



E 
o 
ra 

> -0.05 



-0.1 



X 

\ ^ 

\ ^ 


X / 
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100 
Ag 



FIG. 6. AgTi (Silver - Titanium) ground state convex hull. 
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FIG. 5. AgPd (Silver - Palladium) ground state convex hull. 
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Ag-Y (Silver - Yttrium). The stability of AgY-B2 
and Ag2Y-Cllf, is confirmed [9,10,102]. In the Ag-rich re- 
gion, we do not determine the stability of Ag5iYi4, since 
we do not have the prototype in our library [103,104]. 
Instead, we find AgaY with DOa structure stable with 
an energy ~22nieV/atom below the energy of a two- 
phase region Ag^Ag2Y. In the real system, the presence 
of Ag5iYi4 probably makes structure DOa metastable. 
In the Y-rich region we find AgY2-C37 degenerate with 
the two-phase region Y<->AgY (within ^1.7meV/atom) 
as shown in figure (7). 



Ag-Y system 


Low Temperature Phases comparison chart 


Composition 
% Ag 


Experimental 
(Massalski [9]) 


Ab initio 
result 


33.3 


two-phase region 
above 200° C 


AgY2-C37/tie-line 


50 


B2 


AgY-B2 


66.6 


Clli, 


AgaY-Cllb 


75 


two-phase region 


AggY-DOa 




above 200° C 


(uncertain) 


78.5 


Ag5iGdi4 


unavailable 




20 40 60 80 

Atomic Percent Silver 



100 
Ag 



FIG. 7. AgY (Silver - Yttrium) ground state convex hull. 



Ag-Zr (Silver - Zirconium). The phase diagram of 
Ag-Zr is known partially, and it has been estimated from 
thermodynamic properties [9,10,105,106]. The stability 
of the two known phases AgZr-Bll (7-CuTi prototype) 
and AgZr2-Cll6 is confirmed by our calculation. In ad- 
dition to the known intermetallic compounds, we find a 
stable phase Ag2Zr, which is not present in Massalski [9]: 
the C6 and C32 structures are degenerate. To conclude, 
we find AgZrs-FCCj^^g degenerate with the two-phase 
region Zr^AgZr2 (within ^1.3meV/atom). To address 
the degenerate structures, we further investigate Ag2Zr 
with PAW-GGA potentials, as described in Section (III). 
With PAW, C32 is the most stable compound and C6 is 
2meV/atom above C32. 



Ag-Zr system 


Low Temperature Phases comparison chart 


Composition 
% Ag 


Experimental 
(Massalski [9]) 


Ab initio 
result 


25.0 


two-phase region 
estimated 
above 700° C 


AgZrg 
FCC|J°'^ /tie-line 


33.3 


Cll6 


AgZra-Cllb 


50 


Bll 


AgZr-Bll (7CuTi) 


66.6 


two-phase region 
estimated 
above 700° C 


AgaZr 
C6/C32 (us-lda) 
C32 stable (paw-gga) 

C6^2meV/at. 
above C32 (paw-gga) 




40 60 
Atomic Percent Silver 



FIG. 8. AgZr (Silver - Zirconium) ground state convex hull. 
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Al-Sc (Aluminum - Scandium). The experimental 
phase diagram has compomids at compositions: AI3SC, 
AI2SC, AlSc, and AISC2 [9,10,107-109]. The ab initio 
technique confirms the stability of AI3SC-LI2, AI2SC-CI5, 
AlSc-B2, and AISC2-B82. In the Sc-rich region of the 
phase diagram, we find a new hexagonal stable phase 
AISC3-DO19, which is not present in Massalski [9]. Only 
very limited experimental data is available for this side 
of the phase diagram. 

Other ah initio studies, relevant for this system, can 
be found in references [110-115]. 
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FIG. 9. AISc (Aluminum - Scandium) ground state convex 
huU. 



Au-Cd (Gold - Cadmium). Several compounds 
have been reported for the system Au-Cd, and the 
phase diagram is known with reasonable accuracy 
[9,10,41,117,119,123-126]. We confirm /3'-AuCd-B19 and 
the high-temperature phase /3-AuCd-B2 (with energy 
'--^GmeV/atom higher than B19). At 25% Cd compo- 
sition, we do not find the compound aiAusCd-AgsMg 
because we do not have the proper prototype in our li- 
brary. Instead, we find structures DO24, DO19, ^-^d AI3PU 
(C03V) with degenerate energies. Such structures have 
hexagonal lattices and they are candidates for the phase 
field a2 which has been reported to be a long-period su- 
perstructure with hP lattice [9]. At 33% Cd composi- 
tion, we find three phases Au2Cd-C37, Au2Cd-MoPt2, 
and Au2Cd-C49 with energies 4meV/atom, 14meV/atom 
and 16meV/atom above the tie-line Au3Cd <->AuCd, re- 
spectively. At 50% Cd composition, the prototype of the 
phase /9"-AuCd is not known, but it is reported to have a 
hexagonal structure. For /?", we suggest two candidates: 

AuCd-Llo and AuCd-FCC{^2B2 (CH "40" in reference 
[262]), which have energies higher by 5.4meV/atom and 
5.5meV/atom with respect to B19. None of these struc- 
tures has hexagonal lattice type, and all our hexagonal 
prototypes have much higher energies. 



Au-Cd system 
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Composition 
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Experimental 
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Ab initio 
result 


~25 


ai-AgsMg 


unavailable 
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a2-hP? (unknown) 


AusCd 
DO24/DO19/AI3PU 
hP candidates for 02 


33.3 


two-phase region 
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Au2Cd-C37, MoPt2, C49 
above tie-line 


~47.5 
~50 
43 to 57 


/3'-B19 
P" unknown 
/3-B2 (high T) 


AuCd-B19 

Llo/FCCJ^JM ~5meV/at. 
B2 ~6meV/at. 
above B19 


61.6 to 62.1 


Au3Cd5-D8m 


unavailable 


72 to 76 


e'-unknown 
tentative phase field 


AuCd3-L6o 


~82 


77' -unknown 
tentative phase field 


nothing stable 
at 83.3% 



At concentration ~75% Cd, a stable phase e', with un- 
known structure, is believed to exists [9]. We confirm the 
existence of a stable compound and our best prediction is 
AuCda-LGo. At concentration ~82% Cd another stable 
phase 77' is believed to exist. We are not able to confirm 
77'. In fact at 83.3% Cd, the phase with lowest energy is 
~60meV/atom higher than the two tie-line AuCda^Cd. 
To address the degenerate structures, we further investi- 
gate Au^Cd with PAW-GGA potentials, as described in 
Section (III). Also with PAW, DO24,, Al^Pu and DOig 
remain degenerate. 
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FIG. 10. AuCd (Gold - Cadmium) ground state convex hull. 



Au-Nb (Gold - Niobium). The phase diagram of 
the system Au-Nb is based on three known compounds 
(AuaNb, AuaNba, AuNbg) [9,10,47,117,119,127-141]. 
We do not confirm the stabihty of phase Au2Nb-C32 
[9,10,47]. Instead of C32, we find a stable com- 
pound Au2Nb-C6 and C32 with energy higher by 
^12meV/atom with respect to C6. We can not conclude 
anything about stability of Au2Nb3 since our library does 
not contain any A2B3 prototype. The Nb-rich side of 
the experimental phase diagram is drawn with the as- 
sumption that AuNba-AlS is stable at low temperature 
(--500°C) [9]. Our calculations suggest that AuNbs- 
A15 might not be stable at low-temperature as it is 
'^7meV/atom above the two phase field AuNb2^Nb. In 
addition, we find a stable phase AuNb2 at 66.6% Nb com- 
position, with bcc parent lattice, and space group Fmmm 



[Oil] 
AB2' 



has AB2 

,[011] 



#69. The compound, labeled as BCC 
stacking along [Oil] direction. If the prototype BCC[1^^2 
were not included in the calculation, the phase AuNbs- 
A15 would be stable, as shown in figure (11). To con- 
clude, we think it might be worthwhile to reconsider the 
thermodynamic modeling of the Nb-rich side of the phase 
diagram with the inclusion of AuNb2-BCCj^]j^2: to ob- 
tain the proper temperature and concentration ranges of 
phase A15. To address the disagreements with the exper- 
imental results, we further investigate Au2Nb and AuNb^ 
with PAW-GGA potentials, as described in Section (III). 
For the compound Au2Nb C32 is the most stable com- 
pound and C6 is higher by 3.2meV/atom. For composi- 
tion AuNbs, A 15 is still unstable, being ^6.5meV/atom 
above the tie-line AuNb2^Nb. The disagreement at 
composition AuNbs is further discussed in Section (IV). 
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A15 is ^6.5meV/at. 
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See Section (IV). 
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Au-Pd (Gold - Palladium). Two ordered com- 
pounds have been suggested for the system Au-Pd in the 
fee solid solution (Au,Pd): AugPd and AuPda [142-145]. 
The low-temperature part of the phase diagram has been 
constructed with the addition of a phase field AuPd that 
is believed to exist below ~100°C [9,10,117,119]. In the 
Au-rich part of the phase diagram, we confirm the stabil- 
ity of a compound AusPd, but we find three structures, 
DO23, DO22 and LI2, with degenerate energy. At compo- 
sition Au4Pd, we find the structure Dl^ to be degener- 
ate with the tie-line Au^AuaPd. At composition 33.3% 
Pd, we find a stable compound Au2Pd, but again, two 
structures C49 and C37 are degenerate. At composition 
50%, we confirm the existence of a stable phase AuPd. 
The prototype is not known. Our best guess is AuPd- 
FCC{^2B2 (CH "40" in reference [262]), and LIq with an 
energy higher by 8meV/atom with respect to FCC[^2B2- 
Finally, in the Pd-rich part of the phase diagram, we find 
a compound AuPda to be degenerate with the tie-line of 
the two-phase region AuPd^Pd. For this compound, 
three structures, DO23, DO22 and LI2, have degenerate 
energies. The considerable degeneracy of the ground 
states is indicative of week effective interactions and con- 
sistent with the significant miscibility of the two elements 
inside each other. The near degeneracy of structures 
such as DO23, DO22 and LI2, which are related to each 
other by periodic antiphase boundaries, indicates that 
more complicated Long Period Superstructures (LPS) 
might be present. To address the small energy differ- 
ences between structures, we further investigate Au^Pd, 
Au2Pd, AuPd, and AuPds with PAW-GGA potentials, 
as described in Section (III). For the compound Au^jPd, 
DO23 is the most stable compound and DO22 and LI2 are 
higher by 9.4meV/atom and 12.7meV/atom, respectively. 
For the compound AuPd^, DO23 is the most stable com- 
pound and LI2 and DO22 o,re higher by 1.8meV/atom and 
4--8meV/atom, respectively. Also with PAW, for the com- 
pound Au2Pd, C37 and 0^9 remain degenerate. For the 
compound AuPd, FC(3^g2 is the most stable compound 
and LIq is higher by 10.8meV/atom. 

Another ab initio study, relevant for this system, can 
be found in reference [60]. 
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FIG. 12. AuPd (Gold - Palladium) ground state convex hull. 



Au-Sc (Gold - Scandium). The phase diagram of 
the system Au-Sc is known only in the Au-rich region, 
and a total of three intermetaUic compounds have been 
reported: AU4SC, AuaSc, and AuSc [9,10,119,153-158]. 
At 20% Sc, we confirm Au4Sc-Dla. At 25% Sc, we 
find a two-phase region, in agreement with experiments. 
In fact, AusSc-DOq, the structure with lowest energy, 
is higher by ~7meV/atom with respect to the tie-line 
Au4Sc-Dla^Au2Sc-Cllfc. At 33% and at 50% concen- 
trations of Sc, we find two stable phases Au2Sc and AuSc, 
but in both cases two structures are degenerate Cllb 
and MoPt2 for Au2Sc, and B2 and B19 for AuSc. The 
experimental results indicate Au2Sc-Cllfc and AuSc-B2. 
The AuSc2-C37 (Co2Si prototype) compound has been 
speculated to exist by similarity with other Au-(Rare 
Earth)2 systems [153]. Our ab initio method confirms 
AuSc2-C37. To address the degenerate structures, we fur- 
ther investigate Au2Sc and AuSc with PAW-GGA poten- 
tials, as described in Section (III). For compound Au2Sc, 
with PAW, Au2Sc-Cllb is the most stable compound and 
Au2Sc-MoPt2 is higher by l.OmeV/atom. In addition, 
for compound AuSc, Au2Sc-B2 is the most stable com- 
pound and AU2SC-BI9 is higher by 2.2meV/atom. While 
these energy differences are very small, both calculations 
are in agreement with experiments. 
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Au-Ti (Gold - Titanium). Five ordered compounds 
have been reported for the system Au-Ti at low tem- 
perature: Au4Ti, Au2Ti, aAuTi, /3AuTi, and AuTa 
[9,10,48,101,160-169,311,283,321,298]. We confirm the 
stabihty of aAuTi-Bll, Au4Ti-Dla and AuTis-AlS. At 
25% Ti, we find a two-phase region, in agreement with ex- 
periments. In fact, AuaTi-DOa, the structure with lowest 
energy, is higher by '^llmeV/atom with respect to the 
tie-line Au4Ti^->Au2Ti. For compound Au2Ti, the ex- 
perimental phase is Cllft. At such composition we find 
two structures with degenerate energies: Au2Ti-Cllfc and 
Au2Ti-MoPt2. We do not confirm an off-stoichiometry 
stable phase /3AuTi-B19, because we do not have such 
prototype in our library. The stoichiometric B19 has an 
energy ~50meV/atom above aAuTi-Bll. At 42.8% con- 
centration of Ti, we find a previously unknown compound 
Au4Ti3-Cu4Ti3. This compound is degenerate with the 
tie-line of the two-phase field Au2Ti^AuTi, therefore its 
existence is uncertain. To address the degenerate struc- 
tures Au2Ti-Cllb/MoPt2 and the energy difference be- 
tween aAuTi-Bll and [3AuTi-B19 phases, we further in- 
vestigate Au2Ti and AuTi with PAW-GGA potentials, as 
described in Section (III). With PAW, Au2Ti-Cllb is 
the most stable compound and Au2 Ti-MoPt2 is higher 
by S.lmeY /atom. For phase AuTi, AuTi-Bll is the 
most stable compound (as in the ultrasoft pseudopoten- 
tial case), and AuTi-B19 is higher by 53.6meV/atom. 
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FIG. 14. AuTi (Gold - Titanium) ground state convex hull. 



Au-Y (Gold - Yttrium). Four intermetallic com- 
pounds have been reported for the system Au-Y: AuaY, 
AuaY, AuY, and AuYg [9,10,119,155-158,170-173]. 
With our ab initio method, we confirm stability of AuaY- 
DOa and AU2Y-CII6. At equal composition, AuY, we 
do not confirm B2: we find a stable B33 and phase B2 
higher by '^26meV/atom. No experimental compound is 
reported for composition AuY2 , though we obtain AuY2- 
C37. At 75% composition of Y, we do not find a stable 
AUY3 compound reported experimentally [10,156]: the 
structure with lowest energy is AUY3-DO11 with an en- 
ergy higher by ^30meV/atom with respect to the tie- line. 
We further investigate AuY with PAW-GGA potentials, 
as described in Section (III). With PAW, B33 is the most 
stable compound and B2 is '^25meV/atom higher than 
B33. The disagreement at composition AuY is further 
discussed in Section (IV). 
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FIG. 15. AuY (Gold - Yttrium) ground state convex hull. 
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Au-Zr (Gold - Zirconium). Several ordered com- 
pounds have been reported for the system Au-Zr: Au4Zr, 
AusZr, Au2Zr, AuioZry, Au4Zr5, AuZr2, and AuZra 
[9,10,163,174-177]. Our ab initio method confirms the 
stabihty of AusZr-DOa, Au2Zr-Cllfc and AuZr3-A15. 
At composition AuZr2, Massalski and the Pauhng file 
report Cllb and CuZr2 type compounds, respectively 
[9,10]. Cllb (MoSi2 prototype) and CuZr2 are very sim- 
ilar structures, both tetragonal (tI6) and with the same 
space group #139. Cll;, belongs to the bcc superstruc- 
ture family while CuZr2, a slightly distorted version of 
Cllb, is a closed packed structure [10], therefore, within 
the formalism of the atomic environments, they repre- 
sents two different prototypes [10]. Our library contains 
both prototypes. However, in our calculations, AuZr2- 
Cllb and AuZr2-CuZr2 are both unstable with degener- 
ate energies higher by ~20meV/atom with respect to the 
tie-line AuZr<-^AuZr3. We can not say anything about 
Au4Zr, AuioZrr and Au4Zr5, because our library lacks 
prototypes at the proper concentrations. At AuZr com- 
position we find two structures with degenerate ener- 
gies: Bll and FCC|^2B9- This prediction is unreliable 
because such compounds appear in the two-phase re- 
gion AuioZr7<-^-Au4Zr5, that we are not able to describe 
properly. At 42.8% concentration of Zr, we find Au4Zr3- 
Cu4Ti3. As before, this prediction is unreliable because 
it appears in the two-phase region Au2Zr^AuioZr7, that 
we are not able to describe. To conclude, the inclusion 
of prototypes of AuioZr7 and Au4Zr5 in our calculations, 
might change the stability of the predicted AuZr and 
Au4Zr3 phases. To address the degenerate structures, 
we further investigate AuZr2 with PAW-GGA potentials, 
as described in Section (III). With PAW, AuZr2-CuZr2 
and AuZr2-Cllb are still degenerate, but become sta- 
ble with an energy 6.6meV/atom lower than the tie-line 
AuZr^AuZr^ . 
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FIG. 16. AuZr (Gold - Zirconium) ground state convex hull. 
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Cd-Nb (Cadmium - Niobium). The phase diagram 
for the system Cd-Nb is not known [9]. It has been re- 
ported that CdaNb has LI2 prototype [178]. However, we 
are not able to find any stable compound, so we propose 
that the system is not compound-forming, but will dis- 
play low-temperature immiscibility. In our calculations, 
Cd3Nb-Ll2 is ~70meV/atom higher in energy than the 
phase separation of Cd^Nb. In addition, the structure 
with the lowest formation energy is a bcc superstructure 
at composition CdNb and energy ^^SSmeV/atom higher 
then Cd<->Nb. To address the disagreement with the ex- 
perimental result, we further investigate the relevant com- 
pounds with PAW-GGA potentials, as described in Sec- 
tion (III). With PAW, CdsNb-Lh is >100meV/atom 
higher in energy than the tie-line Cd^Nb. The structure 
with the lowest formation energy is a bcc superstructure 
at composition CdNb^, and energy ^58meV/atom above 
Cd<^Nb. It is unlikely that the ab initio calculations 
have such a large errors. The disagreement at composi- 
tion CdsNb is further discussed in Section (IV). 
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See Section (IV). 



Cd-Pd (Cadmium - Palladium). Four Cd-Pd 
compounds have been identified experimentally at low- 
temperature (7, 71, 7' and (3i [179,180,9]). However, the 
phase boundaries are unknown. We confirm /3i-CdPd- 
LIq. In the Cd-rich part of the phase diagram, we find 
CdaPd (near 7' at 74%Cd). We are not able to determine 
the exact structure of CdaPd but our best candidates are 
DO19, DO24, NbPda-type, and AlgPu-type (C03V). In the 
Pd-rich part of the phase diagram, we find a stable phase 
CdPda. As before, we cannot determine its structure, 
precisely. Our guesses are DO22 and NbPda-type, which 
have degenerate energies. At concentration 33.3% Cd, we 
find a new compound CdPd2-C37 which is not present in 
Massalski [9]. To address the degenerate structures, we 
further investigate CdPd^ and Cd-^Pd with PAW-GGA 
potentials, as described in Section (III). Gompounds 
GdPd'i-D022 and GdPd^-NbPd^ remain degenerate with 
PAW (the energy difference is l.lmeV/atom). In addi- 
tion, for Cd^Pd composition, NbPd^/Al^Pu remain de- 
generate compounds, and DO22, DOig, DO24, are higher 
by 4-7meV/atom, 5.6meV/atom, 6.7meV/atom, respec- 
tively. 
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FIG. 17. CdPd (Cadmium - Palladium) ground state convex 
hull. 
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Cd-Pt (Cadmium - Platinum). The phase dia- 
gram of the system Cd-Pt is partiaUy known [181,9], with 
several compounds of unknown structure. We confirm 
the stability of phase a'j^-CdPt-Llo- In the Cd-rich part 
of the phase diagram, the prototypes of Cd2Pt and 71 
are unknown. For Cd2Pt, our best guess are two degen- 
erate structures, Cd2Pt-C37 and Cd2Pt-C16, which are 
also degenerate with the hull. For 71, at stoichiomctry 
CdaPt, we find two degenerate structures DOn, DOq and 
D022- In the Pt-rich part of the phase diagram, we do 
not find a stable compound a'-CdPt3-Ll2. Instead, we 
find a stable orthorhombic CdPta with fee superstructure 
and Cmmm #65 space group. The prototype, labeled 
as CdPtg'"*", is described in Appendix (XI). In addi- 
tion, CdPt3-Ll2 is found with an energy '--^25meV/atom 
above CdPt^'^"*". To address the degenerate structures 
and the disagreement with experimental results for com- 
pound CdPts, we further investigate CdPt^, Cd2Pt and 
CdjPt with PAW-GGA potentials, as described in Section 
(III). With PAW, CdPts-CdPtl''''*'' is the most stable 
structure and CdPt^-Ll2 higher by ^ lO.^meV/atom. For 
compound Cd2Pt, C16 is the most stable compound and 
C37 is ^33meV/atom higher than C16. For compound 
CdsPt, DOii is the most stable compound and DOa and 
DO22 are higher by r^52meV/atom and ^60meV/atom, 
respectively. The disagreement at composition CdPta is 
further discussed in Section (IV). 
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Cd-Rh (Cadmium - Rhodium). No experimental 
phase diagram is available [9]. With our ab initio tech- 
nique we find two stable compounds, Cd2Rh-C37, and 
CdsRh. Our best guesses for CdsRh are Cd3Rh-D024 
and CdaRh-AlaPu-type (C03V), which have degenerate 
energies (2meV/atom difference). Interesting metastable 
phases are CdsRh-DOig, (11 meV/atom above CdsRh), 
and Cd2Rh-C49 (8meV/atom above CdzRh-CS?). To 
address the degenerate structures, we further investigate 
Cd^Rh with PAW-GGA potentials, as described in Sec- 
tion (III). With PAW, Cd^Rh-Al^Pu is the ground state, 
and Cdj,Rh-D024, has an energy 5meV/atom higher. 
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FIG. 19. CdRh (Cadmium - Rhodium) ground state convex 
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Cd-Ti (Cadmium - Titanium). Little is known 
for the system Cd-Ti at high-temperature [9,10,183,184]. 
At low temperature only two stable intermetallic com- 
pounds have been reported [184]. We confirm the sta- 
bility of Ti2Cd-Cll6. At composition CdTi, Massalski 
and the Pauling file, report Bll and CdTi compounds, 
respectively [9,10]. Bll (7CuTi prototype) and CdTi are 
very similar structures, both tetragonal (tP4) and with 
the same space group #129. Bll belongs to the bcc su- 
perstructure family while CdTi, a slightly distorted ver- 
sion of Bll, is a closed packed structure [10]. Therefore, 
within the formalism of the atomic environments, they 
represents two different prototypes [10]. The ab initio 
calculation is able to relax one structure into the other, 
easily. For CdTi we confirm the experimental results. In 
addition, at concentration 42.8%Cd, we find a compound 
Cd3Ti4-Cu4Ti3, degenerate with the tie- line of the two- 
phase region CdTi2^CdTi. We have not found any other 
stable or metastable compound: we conclude that the low 
temperature part of the phase diagram is complete. Note 
that this system is very similar to Ag-Ti. 



Cd-Ti system 


Low Temperature Phases comparison chart 


Composition 
% Cd 


Experimental 
(Massalski [9]) 


Ab initio 
result 


33.3 




CdTi2-Clli, 


42.8 


two-phase region 
above 200° C 


Cd3Ti4 
Cu4Ti3 /tie-line 


50 


Bll 
CdTi [10] 


CdTi-Bll (CdTi) 



0.02 




HCP 



Ti 



20 40 60 80 

Atomic Percent Cadmium 



100 
Cd 



FIG. 20. CdTi (Cadmium - Titanium) ground state convex 
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Cd-Y (Cadmium - Yttrium). The stability of com- 
pounds YCd-B2, YCd2-C6 and YCdg (CdgEr prototype) 
[9,10,185-187] is confirmed. The Yttrium-rich side of the 
phase diagram is poorly known. In that region we find 
two degenerate phases Y2Cd-C49 and Y2Cd-C37, both 
of which are close to the energy of the two phase region 
Y^CdY (figure (21)). We do not find any other stable or 
mctastable phase. To address the degenerate structures, 
we further investigate Y2Cd with PAW-GGA potentials, 
as described in Section (III). With PAW, C37 is the most 
stable compound and C4-9 is higher by 7.7meV/atom. 

Another ab initio study, relevant for this system, can 
be found in reference [188]. 
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Cd-Zr (Cadmium - Zirconium). The Cd-Zr sys- 
tem is poorly characterized [49,9,189], and compounds 
have been identified at four compositions CdsZr, Cd2Zr, 
CdZr, and CdZr2. Massalski [9] reports CdsZr-Llo- 
AuCu, which we consider a misprint for Cd3Zr-Ll2- 
AuCus, while CRYSTMET and Pauling File databases 
[10,11] report CdsZi-DOa-zSCugTi from reference [190]. 
We confirm the stability of Cd3Zr-Ll2 and CdZr2-CIl;,. 
At composition CdZr, Massalski and the Pauling file, re- 
port Bll and CdTi compounds, respectively [9,10]. Bll 
(7CuTi prototype) and CdTi are very similar structures, 
both tetragonal (tP4) and with the same space group 
#129. Bll belongs to the bcc superstructure family 
while CdTi, a slightly distorted version of Bll, is a closed 
packed structure [10]. Therefore, within the formalism 
of the atomic environments, they represents two differ- 
ent prototypes [10]. The ab initio calculation is able to 
relax one structure into the other, easily. For CdZz, we 
find CdZr-Llo instead of CdZr-Bll, which has an en- 
ergy ~18meV/atom higher than LIq. The prototype 
of Cd2Zr is not known: our best guess is Cd2Zr-Cllf,. 
In addition, we find a stable compound CdZr3-A15 not 
present in Massalski [9], and a metastable LI2 which is 
8.5meV/atom higher than A15. To address the disagree- 
ment with the experimental results for compound CdZr, 
we further investigate LIq and Bll with PAW-GGA po- 
tentials, as described in Section (III). With PAW, CdZr- 
Bll is the most stable compound and CdZr-Llo is higher 
by ISmeV/atom. 
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Mo-Nb (Molybdenum - Niobium). The system 
MoNb has not been studied in great detail and no ex- 
perimental intermetallic compounds have been found 
[9,10,47,197]. A bcc solid solution is reported from 
2400°C up to the melt. We predict four stable com- 
pounds at low temperature: MoNb2-Cllb, MoNb-B2, 
Mo2Nb-Cllb, and MoaNb-DOs, as shown in figure (23). 

Other ab initio studies, relevant for this system, can 
be found in references [90,198-200]. 
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Mo-Pd (Molybdenum - Palladium). The phase 
diagram of this alloy is known from experimental investi- 
gations and thermodynamic calculations [9,10,46,51,201]. 
The only compound is at composition MoPd2 and is 
listed in Massalski as having approximatively the MoPt2 
structure [9,10,201]. Our ab initio method finds that 
the MoPt2 structure is 8meV/atom higher than MoPd2- 
ZrSi2. In addition, we find a stable phase MoPd4-Dla 
(MoNi4 prototype), previously unknown. 
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FIG. 24. MoPd (Molybdenum - Palladium) ground state con- 
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Mo-Pt (Molybdenum - Platinum). The phase di- 
agram of this alloy is known from experimental investiga- 
tions and thermodynamic modeling [9,10,201-210,277]. 
We confirm the stability of MoPt-B19 and MoPt2 
(with MoPt2 prototype). In the Pt-rich part of the 
phase diagram, we find a compound MoPt4-DlQ (MoNia 
prototype), previously unknown. The observed high- 
temperature phase A15, reported as MogPt in Massal- 
ski [9] and as Mo3.2Pto.8 in the Pauling File [10], is 
off-stoichiometry. Our ab initio method finds a stoichio- 
metric MoaPt-AlS to be '^45meV/atom above the two- 
phase region Mo<-^MoPt-B19. For a detailed experimen- 
tal study of such A15 phase, see reference [210]. 
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Mo-Rh (Molybdenum - Rhodium). The phase di- 
agram of the system Mo-Rh is known above ~900°Cand 
it is based on thermodynamic calculations and exper- 
imental results [9,10,51,211-214]. We confirm the sta- 
bility of the two known compounds MoRh-B19 [9] and 
MoRha-CdMgs [10,213]. At concentration 66.6% Rh, we 
find the stable phase MoRh2-C37 (prototype Co2Si), pre- 
viously not identified [9]. 
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Mo-Ru (Molybdenum - Ruthenium). The phase 
diagram of this system is known at medium and high- 
temperature [9,10,215]. There is a cr phase around 38% 
Ru, which is quite common when bcc (Mo) and hep (Ru) 
elements are mixed. We are not able to confirm the a 
phase since we do not have the proper prototype in our 
library. At concentration 75% Ru, we find the stable 
phase M0RU3-DO19 (formation energy ^ 60meV/atom), 
not present in Massalski [9], indicating that it might be 
a low-temperature ordered structure. 

Another ab initio study, relevant for this system, can 
be found in reference [216]. 
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Mo-Ti (Molybdenum - Titanium). To our knowl- 
edge, no intermetallic compounds have been reported 
for the system Mo-Ti [9,10,48,218-220], and it is con- 
sidered a non-compound forming system [10,11]. In 
contrast to the experimental assessment, we find seven 
stable compounds. Four of these have bcc superstruc- 
tures: orthorhombic MoTia, degenerate with the two- 
phase region Ti-f->MoTi2, with space group Immm ^71 
and prototype MoTig*^"*" (Appendix (XI)); orthorhom- 
bic MosTi with space group Immm 7^71 and prototype 
Mo3TiP''°*° (Appendix (XI)); trigonal MoTia-BCCj^^j^', 
and orthorhombic MoTi with space group Imma #74 
and prototype MoTiP''°*° (Appendix (XI)). In addition, 
we find a stable compound Mo2Ti-Cllb, a compound 
Mo4Ti-DlQ degenerate with the tie-line, and finally a 
monoclinic MosTi which has space group C2/m #12 and 
prototype MosTi^™*" (Appendix (XII)). Given the rel- 
atively large value of the formation energy, the experi- 
mental miscibility gap is surprising. This system war- 
rants further study to sort out the apparent discrepancy 
between experiments and ab initio computation. 

Another ab initio study, relevant for this system, can 
be found in reference [221]. 
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Mo-Zr (Molybdenum - Zirconium). Not much is 
known of the system Mo-Zr to produce a reliable phase 
diagram [9,10,49]. We confirm the stability of the only 
known compound Mo2Zr-C15. Experimentally, the Zr- 
rich side of the phase diagram has not been explored 
in detail at low temperature, and it is reported to have 
a two-phase region Mo2Zr-C15*-+(Q!Zr)-A3 above 400°C. 
However, we find two possible new compounds: the or- 
thorhombic phase MoZra with space group Imma #74, 
prototype MoZrg'^"*" (Appendix (XII)), and the mon- 
ocHnic compound MoZrs-MosTi?"'"*" with space group 
C2/m #12, and prototype MogTi^™*" (Appendix (XII)). 
Both compounds MoZrf™*° and MoZrg-MosTiP™*" are 
quasi-degenerate with respect to the two-phase region 
Mo2Zr^(aZr). Hence, the existence of MoZrg'^"*'' and 
MoZrs-MosTi'"""*" is uncertain. It is useful mention- 
ing that we find a metastable orthorhombic phase MoZr 
with space group Imma #74 and prototype similar to 
MoTiP*^"*" (Appendix (XI)) only -2meV/atom above the 
tie-line Mo2Zr<-^(aZr). 

Another ab initio study, relevant for this system, can 
be found in reference [110]. 
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Nb-Pd (Niobium - Palladium). The Pd-Nb phase 
diagram is known with reasonable accuracy in the Pd-rich 
region [9,10,47,259]. The stability of the experimental 
phase NbPd2-MoPt2 and Q;NbPd3-D022 is confirmed. In 
the Nb-rich region, at concentration 33% Nb, we find an 
orthorhombic compound Nb2Pd-BCC[^]j"'^, not present in 
Massalski [9]. The energy of this phase is ~llmeV/atom 
below the tie-line of the two-phase region Nb<-+NbPd2. 

Another ab initio study, relevant for this system, can 
be found in reference [224]. 
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Nb-Pt (Niobium - Platinum). Several intermetal- 
lic compounds have been reported for the system Nb-Pt 
[9,10,47,225-229]. We confirm the experimental phases 
Nb3Pt-A15 and aNbPta-DOa (/JCugTi prototype). The 
prototype of NbPt2 is not reported in Massalski [9] (or- 
thorhombic oI6, with space group Immm). We find 
NbPt2-MoPt2, in agreement with references [10,227,228]. 
At 50% concentration we do not confirm the stability of 
NbPt-B19. Instead of B19 we find NbPt-Llo, and B19 to 
be higher by ~llmeV/atom above LIq. It is possible that 
LIq is a ground state and B19 a high-temperature state. 
We can not say anything about the cr-phase D8f, since we 
do not have such prototype in our library. We further in- 
vestigated NbPt-Llo and NbPt-B19 with PAW-GGA po- 
tentials, as described in Section (III). With PAW, LIq 
is the most stable compound and B19 is ^ 18.5meV/atom 
higher than LIq. 
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Nb-Rh (Niobium - Rhodium). The system Nb- 
Rh is poorly characterized in the range of concentration 
50~80%Rh [9,10,47,226,230,231,51]. We confirm the sta- 
bility of Nb3Rh-Al5, NbRh-Llo. NbRh-B19, which is 
observed as a high-temperature state, is '^27meV/atom 
above Llo, and it is possibly stabilized by entropic ef- 
fects. Note: in Massalski and in the Pauling File [9,10] 
there is no phase at 50% composition, and Nbo.96Rh1.04- 
Llo appears ofF-stoichiometry. At concentration 75% Rh, 
we find the stable phase rj-Al^Pu (C03V) and KNbRha- 
LI2 to be higher by 8meV/atom. Hence, we think that 
ry prevails at low-temperature over k, in contrast with 
the sketched phase diagram of reference [9] . We can not 
say anything about D8fc and ^(Nb2Rh3) since we do not 
have the cr-phase DSfc and any A2B3 prototypes in our 
library. To address the structures with similar energy, 
we further investigate NbRh^ with PAW-GGA potentials, 
as described in Section (III). With PAW, NhRh^-Al^Pu 
is the most stable compound and NbRhz-Ll2 is higher by 
5.3meV/atom. The disagreement at composition NbRhs 
is further discussed in Section (IV). 

Another ab initio study, relevant for this system, can 
be found in reference [232]. 
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Nb-Ru (Niobium - Ruthenium). Very lit- 
tle is known for the alloy Nb-Ru, especially at low- 
temperature [9,10,47,233-237]. We do not confirm the 
stability of any low-temperature compound at compo- 
sition NbRu, but instead find a two-phase field be- 
tween NbsRu and NbRu2, even though a high and 
low-temperature phase of NbRu has been observed 
[233,237]. In the Ru-rich side of the phase diagram, 
we find NbRu3-D024 to be 8meV/atom lower than LI2 
which is suggested experimentally. At 66% Ru we 
find NbRu2-C37 (with oP12-Co2Si prototype). With 
respect to the two-phase field Nb3Ru^NbRu2, the 
structures closest to the tie-line at NbRu composi- 
tion are B19 (~13meV/atom), LIq (~20meV/atom), 
B27 (-23meV/atom), B33 (-39meV/atom), and B2 
('--^45meV/atom). Experiments have not found any other 
stable compound. However, in the Nb-rich side of 
the phase diagram, we find NbaRu-DOa and a mono- 
clinic NbsRu with C2/m #12 space group and proto- 
type Mo5TiP''°*°. The structure NbsRu-MosTi?"'''*" is 
described in Appendix (XII). To address the disagree- 
ments at compositions NbRu and NbRu^ we further in- 
vestigate the relevant compounds with PAW-GGA po- 
tentials, as described in Section (III). For composi- 
tion NbRu, the PAW-GGA result is substantially dif- 
ferent from the US-LDA result. GGA gives LIq only 
4meV/atom above the tie-line (versus 20meV/atom in 
LDA), whereas B19 was difficult to converge numerically, 
and seems to be > lOOmeV/atom above the tie-line. For 
composition NbRu^, also with PAW, DO24 is the most 
stable compound and LI2 is higher by 2.5meV/atom. Be- 
cause this number is so small, extremely large h-points 
sets and high energy cutoff were used to converge it. The 
disagreements at compositions NbRu and NbRua are fur- 
ther discussed in Section (IV). 

Other ab initio studies, relevant for this system, can 
be found in references [238-240]. 
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FIG. 33. NbRu (Niobium - Ruthenium) ground state convex 
liull. 
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Nb-Tc (Niobium - Technetium). The phase di- 
agram of the system Nb-Tc is not known and only 
one intermetalhc compound, NbTca (metallic and su- 
perconductor) has been reported with aMn structure 
[9,10,47,241-243,365]. The compound NbTcs is classified 
as phase Nbo.15Tco.85 in the Pauling File [10]. We do not 
have aMn in our library of prototypes, and, at such sto- 
ichiometry, we do not find any stable compound. In the 
Tc-rich side of the phase diagram, we find a two-phase 
field NbTc^Tc, as shown in figure (34). At 50% con- 
centration, we obtain NbTc-B2, and in the Nb-rich side 
of the diagram, we find Nb2Tc-Cllf, and an orthorhom- 
bic phase NbaTc with space group Immm #71, bcc su- 
perstructure, and prototype NbaTc^*""*" described in Ap- 
pendix (XI) . Trends of Tc alloys are further discussed in 
Section (VIII). 
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FIG. 34. NbTc (Niobium - Technetium) ground state convex 
hull. 



Pd-Pt (Palladium - Platinum). The low temper- 
ature part of the phase diagram is believed to have a 
miscibility gap at a temperature of about 770°C. This 
miscibility gap is predicted on the basis of the difference 
of mehing points between Pd an Pt [9,10,43,260,261]. In- 
stead of such gap, we find three unknown stable com- 
pounds with fee superstructures (Pd and Pt are both 
fee). We find PdPt-Lli, and two orthorhombic phases, 
PdaPt and PdPta, with space group Cmmm #65 and 
prototypes described in Appendix (XI). The prototypes 
are labeled as PdaPt^™*" and PdPtg™*°. The compound 
PdaPt is degenerate with respect the two-phase fields Pd- 
Al^PdPt-Lli, therefore its existence is uncertain. As 
shown in figure (35), all the stable phases have small for- 
mation energy (< 50meV/atom) making them difficult 
to determine experimentally. However, we believe the 
experimental phase diagram to be in error. We further 
investigate PdPt with PAW-GGA potentials, as described 
in Section (III). With PAW, PdPt-Lli is the most stable 
compound and PdPt-Ll^ is higher by 5.5meV/atom. Our 
results are in disagreement with previous FLAPW-LDA 
calculations, where Llo is found to be the most stable 
compound [262]. 

Other ab initio studies, relevant for this system, can 
be found in references [13,62]. 
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Pd-Tc (Palladium - Technetium). The phase dia- 
gram for the system Pd-Tc is sketched based on the sohd 
solubihty data [9,10,43,365,368,275,346]. Experimental 
results report two solid solutions (fee Pd-rich) and (hep 
Tc-rich) with a two-phase region in between. No inter- 
metallic compounds have been reported [9,10]. However, 
we find one stable phases PdTca-DOig. In addition, at 
50% concentration, we find a hep superstructure (trigo- 
nal lattice, hP4, space group P3ml) which has an energy 
'^3meV/atom higher than the tie- line of the two-phase 
region PdTca^Pd. Trends of Tc alloys are further dis- 
cussed in Section (VIII). 
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FIG. 36. PdTc (Palladium - Technetium) ground state con- 
vex hull. 



Pd-Ti (Palladium - Titanium). This system 
has been subject of confiicting results for several years 
[9,10,169,276-291]. Our ab initio method confirms the 
stability of the compounds PdTi2-Cllh, PdTi3-A15, and 
Pd3Ti-D024. Near 80% Pd, Long Period Superstructures 
(LPS) modulations of LI2 are observed [279]. While we 
do not have such off-stoichiometric LPS, we find LI2 to 
be only 6meV/atom above DO24 at PdsTi composition. 
The low energy difference between LI2 and DO24 indi- 
cates that it may be easy to form antiphase boundaries 
and LPS near this composition. We find Pd2Ti-MoPt2 
which is an orthorhombic distortion of Gilt (MoPt2 is 
orthorhombic, while Cllf, is tetragonal). We find Pd2Ti- 
C49 and Pd2Ti-Cllh to be higher by 3meV/atom and 
14meV/atom above MoPt2, respectively. At low tem- 
perature, at concentration 50%, we find a stable com- 
pound aTiPd, but two prototypes are degenerate: LIq 
and B19 (which is reported experimentally). We can not 
find the reported phase PdTi4-Al5 (ofF-stoichiometry) 
[9], which we think should appear inside the two-phase 
region of Ti ~ PdTia, or at composition PdTia. While 
ofi^-stoichiometric compounds are obviously possible, this 
usually goes together with significant width of the single- 
phase field. Hence, we concluded (maybe erroneously) 
that the placement of A15 at composition PdTi4 in ref- 
erence [9] is likely a typographical error. To address the 
degenerate structures, we further investigate a TiPd with 
PAW-GGA potentials, as described in Section (III). With 
PAW, PdTi-B19 is the most stable compound and PdTi- 
LIq is higher by lOmeV/atom. 

Another ab initio study, relevant for this system, can 
be found in reference [292]. 
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FIG. 37. PdTi (Palladium - Titanium) ground state convex 
hull. 



Pd-Y (Palladium - Yttrium). The phase diagram 
for the system Pd-Y is Icnown with reasonable accuracy. 
Several intermetallic compounds have been reported, but 
not all the structures have been determined experimen- 
tally [9,10,293-296]. We confirm the stability of the com- 
pounds Pd3Y-Ll2 and PdYs-DOn. The prototype of 
aPdY is not known [293]. Our best guesses are PdY- 
B27 and PdY-B33 (CrB), which are degenerate in the 
calculation. We also find PdY2-C37, which occurs in 
a concentration between two known compounds Pd2Y3 
and Pd2Y5 and is very close to the tie-line of the two- 
phase field PdYa^PdY (5meV/atom) Hence, if Pd2Y3 
and Pd2Y5 were to be included, C37 would likely not be 
stable. We cannot check this prediction since we do not 
have Pd2Y3 and Pd2Y5 prototypes in the set of calcu- 
lations. To address the degenerate structures, we further 
investigate PdY with PAW-GGA potentials, as described 
in Section (III). With PAW, PdY-B27 is the most stable 
compound and PdY- B 33 is higher by 3.1meV/atom. 
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FIG. 38. PdY (Palladium - Yttrium) ground state convex 
hull. 
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Pd-Zr (Palladium - Zirconium). The experimen- 
tal phase diagram for the system Pd-Zr is based on 
hmited information [9,10,297-299,282,300-304]. We are 
able to confirm stability of Pd3Zr-D024 and PdZr2-Cllb. 
The stable phase of PdZr is reported to be CrB (B33) 
[10,304]. At that composition, we find two degenerate 
structures: PdZr-B27 or PdZr-B33 (CrB). In the Zr- 
rich side of the phase diagram, we find PdZra-FCCj^^^]. 
However, the energy difference between PdZrs-FCCj^^g 
and the two-phase region Zr<-!-PdZr2-Cllf,, is very small 
(~6meV/atom). Hence, the existence of compound 
PdZrs-FCCj^^"';] remains uncertain. At low-temperature, 
instead of a stable phase with stoichiometry Pd2Zr, we 
find the two-phase field PdZr<->Pd3Zr-D024. The least 
unstable phases are Pd2Zr-C49, Pd2Zr-MoPt2, Pd2Zr- 
Cllf), with energies ^18meV/atom, ~24meV/atom, 
~26meV/atom above the tie-line of the two-phase re- 
gion, respectively. To address the degenerate structures 
and the erroneous two phase region PdZr^Pd^Zr, we 
further investigate PdZr and Pd2Zr with PAW-GGA po- 
tentials, as described in Section (III). With PAW, PdZr- 
B33 is the most stable compound and PdZr-B87 is higher 
by 3.2meV/atom. In addition, Pd2Zr-Cllb is degener- 
ate with the tie line PdZr^Pd^Zr, and MoPt2 and C49 
are higher by r^l.SmeV/atom and r^l.lmeV/atom with 
respect to Cllb (also the tie-line PdZr^Pd^Zr has been 
recalculated with PAW-GGA potentials). 
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FIG. 39. PdZr (Palladium - Zirconium) ground state convex 
hull. 



Pd-Zr system 


Low Temperature Phases comparison chart 


Composition 
% Pd 


Experimental 
(Massalski [9]) 


Ab initio 
result 


25 


two-phase 
region 


PdZr3-FCC|^°'^] 
(uncertain) 


33.3 


Cllt 


PdZr2-Clli, 


50 


disorder fee [9] 
CrB(B33-TlI) 
[10,304] 


PdZr- B27/B33 (us-lda) 
B33 stable (paw-gga) 

B27^3.2meV/at. 
above B33 (paw-gga) 


66.6 


Cllfc 


two-phase region (us-lda). 
C49~18meV/at. 
MoPt2~24meV/at. 
Clli,~26meV/at. 
above the tie-line (us-lda). 
Gil}, /tie-line (paw-gga). 
MoPt2^1.3meV/at., 
C49r^7.1meV/at. 
above tie-line (paw-gga). 


75 


DO24 


Pd3Zr-D024 



35 



Pt-Rh (Platinum - Rhodium). The experimen- 
tal phase diagram of the system Pt-Rh is similar the 
system Pd-Pt. The low-temperature part of the phase 
diagram is believed to have a miscibility gap at a tem- 
perature of about 760°C [9,10,43,263,305]. Instead of 
the gap, we find several stable phases, all with fee super- 
structure (Pt and Rh are both fee), similar to previous 
FLAPW-LDA calculations [262]. We find Pt4Rh-Dla, 
Pt3Rh-D022, PtRh2-C49, PtRh3-D022, PtRh4-Dla, and, 
at 50% concentration, PtRh-FCC{^2B2 (CH "40" in refer- 
ence [262]). As shown in figure (40), all the stable phases 
have very small formation energy (< 30meV/atom) indi- 
cating that they may disorder at relatively low tempera- 
ture. 

Another ab initio study, relevant for this system, can 
be found in reference [62]. 
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FIG. 40. PtRh (Platinum - Rhodium) ground state convex 
hull. 



Pt-Ru (Platinum - Ruthenium). Only one 

compound has been found for the system Pt-Ru 
[9,10,306,307]. At low-temperature, the phase dia- 
gram reported in IVIassalski [9] has Platinum-rich and 
Ruthenium-rich solid solution with large solubilities of 
the other element, and a two-phase region for concen- 
tration between ~70% to ~80% of Ruthenium. How- 
ever, recent X-ray diffraction experimental work reported 
the existence of a fee phase at 50% composition, with 
unknown prototype [307]. For PtRu, our prediction is 



PtRu-FCC 



[001] 
A2B2- 



At 25% Ru composition, we find a 

•,[001] 



stable phase PtsRu-FCC^^g, degenerate with the two- 
phase field Pt^PtRu. Hence, the existence of compound 
PtaRu-FCCj^^^l remains uncertain. To our knowledge, 
PtRu is the first known system where the prototype 
structure FCcJ^2B2 would be stable. 
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41. PtRu (Platinum - Ruthenium) ground state convex 
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Pt-Tc (Platinum - Technetium). The phase dia- 
gram for the system Pt-Tc has been determined from ex- 
perimental sohd solubihty data [9,10,43,365,368]. No in- 
termetalhc compounds have been reported [9]. However, 
we find two stable phases: PtsTc-FCCl^^g and PtTca- 
DO19. PtTcs-DOig appears in the composition range of 
a two-phase region Pt-Al and Tc-A3, that is present at 
temperatures higher than ^ 1000°C. Trends of Tc alloys 
are further discussed in Section (VIII). 
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FIG. 42. PtTc (Platinum - Technetium) ground state convex 
hull. 



Pt-Ti (Platinum - Titanium). Not much is 
known of Pt-Ti system to produce a precise phase di- 
agram [9,10,48,285,290,308-311,283,312]. Some inter- 
metallic compounds are reported [308-311,283,312,288]. 
We confirm the stabihty of phases PtTi3-A15, aPtTi- 
B19, Pt3Ti-D024. B19 is stable at 50% composition, 
and LIq and B33 are higher by ^20meV/atom and 
~30meV/atom, respectively. At concentration 25%Ti, 
we confirm the existence of 7-LI2, which has an energy 
that is ^5meV/atom higher than Pt3Ti-D024. However, 
LI2 has been reported to be stable away from stoichioni- 
etry < 25% Ti [9]. At composition 33% Ti, PtsTi, Mas- 
salski reports a two-phase region above 600° C [9,10]. In- 
stead of the two-phase field, we find a stable compound 
Pt2Ti: two structures, C49 and C37, are degenerate. We 
cannot say anything about PtgTi-Dla, because our li- 
brary does not contain off-stoichiometry Dl^. To address 
the degenerate structures, we further investigate Ptz Ti 
with PAW-GGA potentials, as described in Section (III). 
With PAW, Pt2Ti-C4-9 is the most stable compound and 
Pt2Ti-C37 is higher by 1.9meV/atom. 

Another ab initio study, relevant for this system, can 
be found in reference [292]. 
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FIG. 43. PtTi (Platinum - Titanium) ground state convex 
hull. 
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Pt-Y (Platinum - Yttrium). The experimental 
phase diagram of the system Pt-Y has been sketched 
by analogy with other Rare Earth-Platinum diagrams 
[9]. Several intermetallic compounds have been reported 
[9,10,43,156,277,285,313-320,324]. The stability of the 
compounds Pt3Y-Ll2, Pt2Y-C15, and PtYa-DOn is con- 
firmed. At 50% concentration, we do not find any sta- 
ble PtY-B27. Instead of B27, we find PtY-B33 (CrB 
prototype) to have lowest energy with B2 and B27 hav- 
ing energies ^50meV/atom above B33. At 66.6% Y 
concentration, we find PtY2-C37 which is the proto- 
type of Co2Si and Ni2Si. Hence our calculations con- 
firm the correct experimental structure. At concentra- 
tion 62.5% Y, we do not find PtaYs-DSg, but the two- 
phase field PtY^PtY2. DSs and WsSia are higher by 
~24meV/atom and ~80meV/atom with respect to the 
tie-line PtY^PtY2 (this problem is solved with PAW- 
GGA potentials). As shown in figure (44), at concen- 
tration 16.6% Y, we find a two-phase region instead 
of the reported stable compound PtsY with unknown 
structure [9]. Our best guess is Pt5Y-D2£i which is the 
least-metastable structure we have at such composition 
(~18meV/atom above the tie-line Pt<-^Pt3Y). We con- 
clude that further experimental and theoretical investi- 
gations are necessary to determine the behavior of PtY. 
To address the disagreement with the experimental re- 
sults for compounds PtY and Pt^ Y^, we further investi- 
gate the relevant structures with PAW-GGA potentials, 
as described in Section (III). With PAW, PtY-B33 is 
stable and PtY-B27 is higher by 60meV/atom. For com- 
pound Pt^ Y^, D8s is stable with an energy 6.9meV/atom 
lower than the tie-line PtY^PtY2, which has also been 
recalculated with PAW-GGA potentials. The disagree- 
ment at composition PtY is further discussed in Section 
(IV). 
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FIG. 44. PtY (Platinum - Yttrium) ground state convex hull. 
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Pt-Zr (Platinum - Zirconium). The system Pt- 
Zr is quite interesting. Our ab initio metlrod con- 
firms the stability of aPtZr-B33 (CrB prototype). Two 
crystal structures have been reported for PtaZr: DO24 
and LI2 [9,10,321,323,322,324,285,302]. We confirm the 
stability of Pt3Zr-D024 and we find LI2 to be higher 
by lOmeV/atom with respect to DO24. In the Zr- 
rich part of the phase diagram, we find two stable 
phases PtZr2-C16 and PtZr3-A15. At concentration 
62.5% Zr, we do not find PtgZrs-DSg, but the two- 
phase field PtZr^PtZr2. WsSia and DSg are higher 
by '--^26meV/atom and ~36meV/atom with respect to 
the tie-line PtZr^PtZr2. To address the phase insta- 
bility at composition Pt^Zri^, we further investigate the 
relevant compounds with PAW-GGA potentials, as de- 
scribed in Section (III). With PAW, at concentration 
Pt^Zr^, there is a two-phase region PtZr'^PtZr2. In 
addition, Pt^Zr^-W^Sis and Pt^Zrf,-D8^ have energies 
higher by 23meV/atom and 26meV/ atoms with respect 
to the tie-line PtZr'^PtZr2, respectively. The disagree- 
ment at composition PtsZrs is further discussed in Sec- 
tion (IV). 
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Rh-Ru (Rhodium - Ruthenium). The phase dia- 
gram for the system Rh-Ru is based on the sohd solubil- 
ity data [9,10,43,326]. Experimental results report two 
solid solutions (fee Rh-rich) and (hep Ru-rich) with a 
two-phase region in between (from 34.5% to 40% atomic 
percent Ruthenium). No intermetallic compounds have 
been reported [325-327] and the system is considered to 
be no compound forming [10]. However, we find two sta- 
ble phases: an orthorhombic oC12 RhRu2P''°*'', with hep 
superstructure and Cmcm #63 space group, and a trig- 
onal hP4 RhRu^"""*", with hep superstructure and P3ml 
#164 space group. Both prototypes are described in Ap- 
pendix (XI). As shown in figure (46), all the stable phases 
have small formation energy (< lOmeV/atom) making 
them difficult to determine experimentally. To better de- 
scribe the stability of this system, we further investigate 
all the structures with negative formation energies with 
PAW-GGA potentials, as described in Section (III). Also 
with PAW, RhRu2P''°*° and RhRuP"^"*" are the most sta- 
ble structures, with formation energies of -8.3meV/atom 
and -8.8meV/atom, respectively. 
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FIG. 46. RhRu (Rhodium - Ruthenium) ground state convex 
hull. 



Rh-Tc (Rhodium - Technetium). The phase di- 
agram for the system Rh-Tc is based on the solid sol- 
ubility data [9,10,43,366,368]. Experimental results re- 
port two solid solutions (fee Rh-rich) and (hep Tc-rich) 
with a two-phase region in between. No intermetallic 
compounds have been reported [366,368]. However, we 
find three stable phases: RhTcs-DOig, RhTc-B19, and 
Ru2Tc-ZrSi2. Trends of Tc alloys arc further discussed 
in Section (VIII). 
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Rh-Ti (Rhodium - Titanium). There are qual- 
itative disagreements about the phase diagram of the 
Rh-Ti system [9,10,48,322,328-332]. We confirm the 
stable phases that where found by all investigators 
[328-331]: RhTi2, aRhTi-Llo, and Rh3Ti-Ll2. Refer- 
ences [10,322,332] report RhTi2-CuZr2 instead of RhTi2- 
Cllfe. The prototype CuZr2 is a distortion of Cllfc and 
has the same lattice type (tetragonal, tI6) and space 
group (I4/mmm #139) [10]. We find RhTis-Cllf, and 
RhTi2-CuZr2 to have degenerate energy. In the Rh-rich 
part of the phase diagram we find a stable phase Rh2Ti- 
C37. However, we do not have RhsTia in our library, 
so C37 might likely be unstable with respect to the two- 
phase field Rh5Ti3<->Ll2. At concentration « 84% Rh, 
we do not find any stable compound, in agreement with 
[330] and in contrast with [328]. To address the degener- 
ate structures RhTi2-Cllb/CuZr2, we further investigate 
RhTi2 with PAW-GGA potentials, as described in Section 
(III). Also with PAW, Cllb and CuZr2 remain degener- 
ate. 

Other ab initio studies, relevant for this system, can 
be found in references [232,292]. 
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Rh-Y (Rhodium - Yttrium). Several com- 
pounds have been reported for the system Rh-Y at low- 
temperature [9,10,43,333,277,334,335,337]. The stabil- 
ity of Rh2Y-C15, RhY-B2, RhYg-DOn is confirmed. At 
composition RhsY, where a D2d structure has been seen 
at high-temperature (but not stable at low-temperature) , 
we find the D2d to be the lowest energy structure (of 
all the structures at that composition), even though it 
is metastable with respect to the phase separation into 
Rh2Y-C15^Rh. At composition RhaY one compound 
has been reported to be stable with prototype CeNia 
and space group PGa/mmc #194 [10,336,338]. We do not 
have such prototype in our library and we do not find any 
stable phase: RhaY-DOig is the least metastable proto- 
type we obtain. RhaY-DOig is higher by 130 meV/atom 
with respect the tie-line, which is at least one order of 
magnitude bigger than the accuracy of the calculations. 
We also find RhY2-C37, which appears in a concentration 
between two known compounds, RhaYs and RhaYy, that 
are not present in our set of calculations. Therefore C37 
might be unstable with respect to the two-phase region 
RhaYs^RhgYy. 
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FIG. 49. RhY (Rhodium - Yttrium) ground state convex hull. 
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Rh-Zr (Rhodium - Zirconium). Although the sys- 
tem Zr-Rh is well known for its superconducting phases 
[9,10,339,342], further investigations are needed to clar- 
ify the stability and presence of intermediate phases 
[340,341]. Our ab initio method confirms the stability 
of RhZr2-C16 and Rh3Zr-Ll2. Massalski does not re- 
port the prototype of the low-temperature phase aRhZr 
[9]. References [342-344] report Q;RhZr-B27. We con- 
firm the stability of aRhZr-B27 (BFe prototype) with no 
other metastable compounds with similar energy. In ad- 
dition, we find three new phases Rh2Zr-C37, KhZi^-Dla 
and RhZrg-FCCl^si' ^^ich are degenerate with the two- 
phase fields RhZr^RhgZr, Zr ^RhZr, and Zr ^RhZr, 
respectively. 

Other ab initio studies, relevant for this system, can 
be found in references [232,345]. 
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FIG. 50. RhZr (Rhodium - Zirconium) ground state convex 
hull. 



Ru-Tc (Ruthenium - Technetium). The phase di- 
agram for the system Ru-Tc is considered to have a con- 
tinuous disordered hep solid solution at low-temperature 
[9,10,365,368,346]. We find three stable ordered phases 
RU3TC-DO19, RuTc-B19, and RuTcg-DOig, all of which 
are hep superstructures. Hence, it is possible that they 
are low-temperature phases of the system Ru-Tc. Trends 
of Tc alloys are further discussed in Section (VIII). 
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FIG. 51. RuTc (Ruthenium - Technetium) ground state con- 
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Ru-Ti (Ruthenium - Titanium). The phase di- 
agram of Ru-Ti is well determinated by several inves- 
tigators. A single low-temperature compound RuTi-B2 
has been reported [9,10,48,285,347,348,350,349,351-354], 
and our ab initio method confirms its stability. In ad- 
dition, we find two stable phases: RuTi2-C49, and or- 
thorhombic RuTis with space group Immm ^^71, bcc su- 
perstructure, and prototype RuTig*^"*" described in Ap- 
pendix (XI). Such compounds are close to the tie- line 
Ti^RuTi. In fact, for C49 and RuTif the formation 
energies are lower by '^34meV/atom and ~37meV/atom 
with respect to the two-phase field Ti<-->RuTi. 
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FIG. 52. RuTi (Ruthenium - Titanium) ground state convex 
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Ru-Y (Ruthenium - Yttrium). Several com- 
pounds have been reported for the system Ru-Y 
[9,10,43,333,355-359]. We confirm the stability of RuYg- 
DOii and Ru2Y-C14. Also, we find RUY2-CI6, which ap- 
pears in a concentration between two known compounds 
RU25Y44 and RU2Y5 that are not in our library of calcu- 
lations. Hence, the existence of RUY2-CI6 is uncertain. 
At concentration 50%, we do not find any stable RuY 
compound, in agreement with [358]. 



Ru-Y system 


Low Temperature Phases comparison chart 


Composition 
% Ru 


Experimental 
(Massalski [9]) 


Ab initio 
result 


25 


DOii 


RuYs-DOii 


28.6 


Ru2Y5-C2Mn5 
mS28 C12/cl 
[10,355,356,358] 


unavailable 


33.3 


two-phase region 


RUY2-CI6 




above 0°C 


(uncertain) 


36.2 


RU25Y44 (unknown) 
oP276 Pnna 
[10,355,359] 


unavailable 


40 


Ru2Y3-Er3Ru2 
hPlO P63/m [360] 


unavailable 


66.6 


RU2Y-CI4 


RU2Y-CI4 




20 40 60 80 100 

Y Atomic Percent Rutlienium Ru 

FIG. 53. RuY (Ruthenium - Yttrium) ground state convex 
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Ru-Zr (Ruthenium - Zirconium). The phase 
diagram of RuZr is known accurately [9,10,362,361]. 
Our ab initio method confirms the stabiUty of the low- 
temperature phase RuZr-B2. In agreement with experi- 
ments, we find no ground state at composition Ru2Zr, 
though the lowest energy structure at that composi- 
tion in our calculations is C14 which appears in the 
phase diagram at high-temperature. Ru2Zr-C14 is higher 
by ~60meV/atom with respect to the two-phase field 
Ru<->RuZr. In the Zr-rich region we find one stable 
compound, HuZr^-Dla, previously unknown. In addi- 
tion, we find three metastable phases: RuZr5, RuZrs, 
and RuZr2-C49, with energies higher by 21meV/atom, 
16.6meV/atom, and 14meV/atom, with respect to 
the two-phase fields RuZr4<->Zr, RuZr<->RuZr4, and 
RuZr^RuZr4, respectively. The structure of RuZrs is 
similar to MoZrl*^"*" (Appendix (XII)), while the struc- 
ture of RuZra is similar to RuTig"^"*" (Appendix (XI)). 

Other ab initio studies, relevant for this system, can 
be found in references [238-240,363]. 
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FIG. 54. RuZr (Ruthenium - Zirconium) ground state convex 
hull. 



Tc-Ti (Technetium - Titanium). The phase dia- 
gram of the system TcTi has been constructed by anal- 
ogy with chemically related systems [9,10,48,364,365]. 
Two intermetallic compound, TcTi-B2 and x ^re re- 
ported [364,365]. We confirm the stability of TcTi-B2, 
but we can not say anything about x since we do not 
have prototypes at composition 85% Ru. In addition, 
we find Tc2Ti-Cllb, TcTi2-C49, and an orthorhombic 
phase TcTig with space group Immm #71 and proto- 
type TcTig''"*" described in Appendix (XI). These inter- 
metallics have large negative formation energy, therefore 
they are expected to be very stable. Trends of Tc alloys 
are further discussed in Section (VIII). 
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FIG. 55. TcTi (Technetium - Titanium) ground state convex 
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Tc-Y (Technetium - Yttrium). Not enough in- 
formation exists in order to construct a phase dia- 
gram for the system Tc-Y [9,10]. Only one intermetal- 
hc compound has been reported: TC2Y-CI4 (Friauf- 
Laves/Frank-Kasper phase) [367,369]. We confirm the 
stabihty of TC2Y-CI4. In addition, we find another sta- 
ble phase TCY3-DO11. Trends of Tc alloys are further 
discussed in Section (VIII). 
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ther discussed in Section (VIII). 
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Tc-Zr system 
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FIG. 57. TcZr (Technetium - Zirconium) ground state convex 
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VIII. TREND FOR TECHNETIUM ALLOYS 

In our set of calculations, we have noticed that the 
phase DO19 appears in systems IVETca where ]V[ is a tran- 
sition metal in the columns on the right of Tc (Tc is 
in column 7B) while DO19 is not present if M is in the 
columns on the left of Tc: DO19 is stable in PdTca, PtTcs, 
RhTc3, RuTc3, and unstable in NbTc, TcTi, TcY, and 
TcZr. 
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X. TABLES OF RESULTS 

A. Experimental compounds in agreement with ab initio solutions 
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»7(Nbi3Rli27)-Al3Pu 


ATI 1 A 1 

NbRli3-Al3Pu 


Ag-Zr 


AgZr2-Clli, 


AgZr2-Cll[, 


Pd-Ti 


PdTi3-A15 at 20% 


PdTi3-A15 at 25% 


Ag-Zr 


AgZr-Bll 


A rv "n -1 -I / rT~>"\ 

AgZr-Bll (7CUT1) 


Pd-Ti 


PdTi2-Cll6 


PdTi2-Clli, 


Al-Sc 


AISC2-B82 


AISC2-B82 


Pd-Ti 


a(TiPd)-B19 


PdTi-B19 


Al-Sc 


AlSc-B2 


AlSc-B2 


Pd-Ti 


Pd3Ti-D024 


Pd3Ti-D024 


Al-Sc 


AI2SC-CI5 


AI2SC-CI5 


Pd-Y 


PdY3-D0ii 


PdY3-D0ii 


Al-Sc 


AI3SC-LI2 


AI3SC-LI2 


Pd-Y 


Pd3Y-Ll2 


Pd3Y-Ll2 


Au-Cd 


p AuCd-B19 


AuCd-B19 


Pd-Zr 


PdZr2-Clli, 


PdZr2-Cll(, 


Au-Cd 


pAuCd-B2 (nigh i) 


AuCd-B2 


Pd-Zr 


FdZr-Bd3 


FdZr-B33 


Au-Nb 


Au2Nb-C32 


Au2Nb-C32 (paw-gga) 


Pd-Zr 


Pd2Zr-Cllj, 


Pd2Zr-Clli,/tie-lme (paw-gga) 


Au-Sc 


Au4Sc-Dl„ 


Au4Sc-Dla 


Pd-Zr 


Pd3Zr-D024 


Pd3Zr-D024 


Au-Sc 


Au2Sc-Clli, 


A -I -I /TV J T^j 

Au2Sc-Clli,/MoPt2 


Pt-Ti 


Pt3Ti-D024 


Pt3Ti-D024 


Au-Sc 


AuSc-B2 


Aubc-B2/B19 


Pt-Ti 


«PtTi-B19 


PtTi-B19 


Au-Ti 


Au4Ti-Dla 


Au4Ti-Dla 


Pt-Ti 


PtTi3-A15 


PtTi3-A15 


Au-Ti 


Au2Ti-Clli, 


A rxT -1 1 /TV T nj. 

AU2 ii-Clli,/MoFt2 


Pt-Y 


Pt3Y5-D88 


l~>j_ "VT" T~VO / ^ 

Pt3Y5-D88 (paw-gga) 


Au-Ti 


aAuTi-Bll 


AuTi-Bll 


Pt-Y 


Pt3Y-Ll2 


Pt3Y-Ll2 


Au-Ti 


AuTi3-A15 


AuTi3-A15 


Pt-Y 


Pt2Y-C15 


Pt2Y-C15 


Au-Y 


AusY-DOo 


AusY-DOa 


Pt-Y 


PtY2-C37 


PtY2-C37 


Au-Y 


Au2Y-Clli, 


Au2Y-Clli, 


Pt-Y 


PtYa-DOii 


PtYs-DGii 


Au-Zr 


AuaZr-DOa 


AusZr-DOo 


Pt-Zr 


Pt3Zr-D024 


Pt3Zr-D024 


Au-Zr 


AugZr-Cllj, 


Au2 Zr-Cllf, 


Pt-Zr 


PtZr-B33 


PtZr-B33 


Au-Zr 


AuZr2-CuZr2/Cll6 


AuZr2-CuZr2 (paw-gga) 


Rh-Ti 


RhTi2-Clli, [9]/CuZr2 [10] 


RhTi2-Cll6/CuZr2 


Au-Zr 


AuZr3-A15 


AuZr3-A15 


Rh-Ti 


oRhTi-Llo 


RhTi-Llo 


Cd-Pd 


CdPd-/3i-Ll(, 


CdPd-Ll() 


Rh-Ti 


Rli3Ti-Ll2 


Rli3Ti-Ll2 


Cd-Pt 


CdPt-a'^-Llo 


CdPt-Llo 


Rh-Y 


RI1Y3-DO11 


RI1Y3-DO11 


Cd-Ti 


CdTi2-Cll6 


CdTi2-Cll6 


Rh-Y 


RhY-B2 


RhY-B2 


Cd-Ti 


CdTi-Bll 


CdTi-Bll (7CuTi) 


Rh-Y 


RI12Y-CI5 


RI12Y-CI5 


Cd-Y 


CdY-B2 


CdY-B2 


Rli-Zr 


RhZr2-C16 


RhZr2-C16 


Cd-Y 


Cd2Y-C6 


C(l2Y-r6 


Rli-Zr 


nRliZr-B27 


RhZr-B27 


Cd-Y 


(:d3Y-Cd3Er 


CkkiY-CkhEr 


Rh-Zr 


Hli3Zr-Ll2 


Rli3Zr-Ll2 


Cd-Zr 


CdZra-Cll;, 


CdZr2-Cll6 


Ru-Ti 


RuTi-B2 


RuTi-B2 


Cd-Zr 


CdZr-Bll 


CdZr-Bll (paw-gga) 


Ru-Y 


RuYs-DOii 


RuYs-DOii 


Cd-Zr 


Cd3Zr-Ll2 


Cd3Zr-Ll2 


R,u-Y 


RU2Y-CI4 


RU2Y-CI4 


Mo-Pt 


MoPt-B19 


MoPt-B19 


Ru-Zr 


RuZr-B2 


RuZr-B2 


Mo-Pt 


MoPt2 


MoPt2 


Tc-Ti 


TcTi-B2 


TcTi-B2 


Mo-Rh 


MoRh-B19 


MoRh-B19 


Tc-Y 


TC2Y-CI4 


TC2Y-CI4 


Mo-Zr 


Mo2Zr-C15 


Mo2Zr-C15 


Tc-Zr 


Tc2Zr-C14 


Tc2Zr-C14 


Nb-Pd 


NbPd2-MoPt2 


NbPd2-MoPt2 









TABLE 5: Experimental compounds in agreement with ab initio solutions (89 entries). 
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B. Experimentally unknown, non-identified or speculated compounds and ab initio predictions 



Experimental unknown/speculated compound ab initio compound 


System 


Experimental result 


Ab initio result 


Ag-Cd 


/3'-bcc ordered 


AgCd-B2/B19/B27 


Ag-Mg 


cF* (unknown) 


AgMg3-D0i9 (hP8) / DOa (oP8) 


Au-Cd 


Au3Cd-a2-liP? (unknown) 


Au3Cd-D024/D0i9/Al3Pu 


Au-Cd 


a" ~ AuCd (unknown) 


AuCd-Llo/FCC^^°!|„ 


Au-Cd 


e'~AuCd3 (unknown) 


AuCd3-L6o 


Au-Pd 


Au3Pd-Ll2 (speculated) 
(wide concentration range) 


Au3Pd-D023 /DO22 /LI2 
Au4Pd-Dla /tic-line 


Au-Pd 


AuPd (speculated and unknown) 


AuPd-FCC''"'' 


Au-Pd 


AuPd3-Ll2 (speculated) 


AuPd3-D023 /DO22 /LI2 /tie-line 


Au-Sc 


AuSc2-C37 (speculated) 


AuSc2-C37 


Cd-Pd 


7'-(unknown) 


Cd3Pd-D0i9/D024/NbPd3/Al3Pu 


Cd-Pt 


7i - (unknown) 


Cd3Pt-DOii/DOa/D022 


Cd-Pt 


Cd2Pt-(unknown) 


Cd2 Pt- C37/C16/tie-line 


Cd-Zr 


cubic (unknown) 


Cd2Zr-Cllf, 


Mo-Pd 


MoPd2 ~MoPt2 


MoPd2-ZrSi2 


Mo-Rli 


M()Rh.■^ (unknown) 


MoRli:i-C(lMK:i 


iNb-1'l 


Nlj Lwu-i)lia,sc ii;giuni — ,£lH.i+,i;-JJ19 


Nbl'l-Li() 


Pd-Ti 


orthorhombic distortion of Pd2Ti-Cll(, 


Pd2Ti-MoPt2 (distortion of Cllt) 


Pd-Y 


aPdY (unknown) 


PdY-B27 


Pt-Ru 


PtRu-FCC (unknown) 


PtRu-FCC 


Pt-Y 


PtsY (unknown) 


Pt5Y-D2d metastable 


Tc-Zr 


TcZr (unknown) 


TcZr-B2 



TABLE 6: Experimentally unknown, non-identified or speculated compounds and ab initio predictions (21 entries). 
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C. Experimental solid solutions, two-phases and "not studied" regions and possible ab initio predictions 



Experimental non-compound <^ ab initio compound 


o y o Li^iii 


J— JA-^Cl lllldlliOjl ICvOLilLi 


blVvVbxJ ICoUlli 


Ag-Au 




r\.^/^r\.\X i^a 1 tit; llllt; 


Ag-Au 


short-range order >950°C 


Ag3 AU-LI2/DO23 /AlsPu/NbPda /DO22 /DO24 


Ag-Au 


short-range order ^950° C 


Ag2Au-C37/MoPt2 


Ag-Au 


short-range order >950°C 


AgAu-Llo 


Ag-Au 


short-range order >950°C 


AgAu2-C37/MoPt2 


Ag-Au 


short-range order >950°C 


AgAu3-Ll2/D022/D023 


Ag-Cd 


fee solid solution 


Ag3Cd-D022/D024 


Ag-Cd 


fee solid solution 


Ag2Cd-C37 


Ag-Cd 


none 


AgCd2-ZrSi2 


Ag-Cd 


hep solid solution 


AgCd3-D0i9 


Ag-Pd 


solid solution >900°C 


AgPd-Lli 


Ag-Pd 


solid solution >900°C 


Ag2Pd-C37 


Ag-Pd 


solid solution >900°C 


AgaPd- LI2/DO22 


Ag-Y 


two-phase region above 200° C 


AgY2-C37/tie-line (uneertain) 


Ag-Y 


two-phase region above 200° C 


AgaY-DOa (uncertain) 


Ag-Zr 


two-phase region above 700° C 


AgZrs-FCC /tie-line (uncertain) 


Ag-Zr 


two-phase region above 700° C 


Ag2Zr-C32 


Al-Sc 


two- phase region above 0°C 


AlSe3-D0i9 


Au-Nb 


two-phase region (ealeulated) 


AuNb2-BCC'?R' 


Au-Pd 


solid solution 


Au2Pd-C49/C37 


Au-Ti 


two- phase region above 500° C 


Au4Ti3-Cu4Ti3/tie 


Au-Y 


not studied 


AuY2-C37 


Au-Zr 


two-pliasc region 


Ati4 Zr3-Cu4Ti3 (unreliable) 


Au-Zr 


L \ V LJIItloV.. 1 v.. IV.Jii 


) 1 1 /KOC'""'^- fiiiivcliahlct 


Cd-Pd 


PH nhase above 100°C 


CdPdQ-DOoo /NhPd-j 


Cd-Pd 


two-phase region above 100° C 


CdPd2-C37 


Cd-Rh 


not studied 


Cd2Rh-C37 


Cd-Rh 


not studied 


Cd3Rh-Al3Pu 


Cd-Ti 


two- phase region above 200° C 


Cu4Ti3/tie-line 


Cd-Y 


not studied/two-phase region 


CdY2-C37 


Cd-Zr 


not studied/two- phase region 


CdZr3-A15 


Mo-Nb 


not studied 


MoNb2-Cllf, 


Mo-Nb 


not studied 


MoNb-B2 


Mo-Nb 


not studied 


Mo2Nb-Cll6 


Mo-Nb 


not studied 


MosNb -DO3 


Mo-Pd 


disorder fee Pd-Al 


MoPd4-Dla 


Mo-Pt 


two-phase region MoPt2<->Pt 


MoPt4-Dla 


Mo-Rh 


two-phase region above 900° C 


Mo2Rh-C37 


Mo-Ru 


disorder hep Ru-A3 above 800° C 


M0RU3-DO19 


Mo-Ti 


not studied/two-phase region above 400°C 


MoTi3'^°'^°/tie-line (uncertain), Appendix (XI) 


Mo-Ti 


not studied/two-phase region above 400°C 


MoTi2-BCCi^gj 


Mo-Ti 


not studied/two-phase region above 400° C 


MoTiP''°«°, Appendix (XI) 


Mo-Ti 


not studied/two-phase region above 400° C 


Mo2Ti-Cllj, 


Mo-Ti 


not studied/two-phase region above 400° C 


Mo3TiJ"'°*°, Appendix (XI) 


Mo-Ti 


not studied/two-phase region above 400°C 


Mo4Ti-Dla/tie-line 


Mo-Ti 


not studied/two-phase region above 400°C 


MosTiJ"""*", Appendix (XII) 


Mo-Zr 


two- phase region above 400° C 


MoZr?'^°*°/tie-line, Appendix (XII) (uncertain) 


Mo-Zr 


two-phase region above 400° C 


MoZr3'^°*°/tie-line, Appendix (XII) (uncertain) 


Mo-Zr 


two-phase region above 400° C 


~ MoTiJ"'°«° (uncertain) 
(~2meV above two-phase reg.) 



TABLE 7a: Experimental solid solutions, two-phases and "not studied" regions and possible ab initio predictions. 
The table continues in the next page. 
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Experimental non-compound <^ ab initio compound 


System 


Experimental result 


Ab initio result 


Nb-Pd 


two-phase region above 700° C 


A1 1 1 

NbaPd-BCC'"^^ 


Nb-Ru 


disorder Nb-A2 


NbsRu Appendix (XII) 


Nb-Ru 


disorder Nb-A2 


NbsRu-DOs 


Nb-Rii 


t\v()-]ilias(~ K^gioii ab(n'(~ 700"C 


NbRii2-C37 


Nb-Tc 


not bLudicxl/ UliivliOWll 


>JI.);iTc'"""', Appendix 1_X1) 


Nb-Tc 


not studied/unknown 


Nb2Tc-Cll6 


Nb-Tc 


not studied/unknown 


NbTc-B2 


Pd-Pt 


two-phase region Pd^Pt 


Pd3Ptf"'°*°/tie-line (uncertain), Appendix (XI) 


Pd-Pt 


two-phase region Pd^ — ^Pt 


PHPf T 1 1 
r lAr L-ij 1 1 


Pd-Pt 


two-phase region Pd^^Pt 


PdPt?J''°", Appendix (XI) 


Pd-Tc 


two-phase region above ~1000°C 


PdTc3-D0i9 


Pd-Y 


two-phase region 


PdY2-C37 (uncertain) 


Pu-Zr 


two-phase region 


PdZirs-l" (^O^^^^ (uncertain j 


r^t-rtn 


two-pnase region rt<-^ixn [yj, Uia [^DzJ 


0+ T3U T-\1 


Jr^t-Kn 


two-phase region Pt^Rh 


Pt3Kh-L)022 


rt-txil 


two-phase region Pt^Rh 


FtKn-lN br^ 


rt-txil 


two-phase region Pt^Rh 


PtKh2-04y 


r^t--tln 


two-phase region Pt-*— ^Rh 


r>+r>u T^n 
FtKn3-UU22 


r^t-rln 


two-pnase region Pt*->Kn [yj, JJia [^bzj 


TD+OT^ T— 11 
Ftrv,n4-Uia 


Pt-Ru 


disorder Pt-Al 


Pt3Ru-FCC'Jp':| /tie- line (uncertain) 


Pt-Tc 


disorder Pt-Al 


Pt3Tc-FCC'°°^^ 


Pt-Tc 


two- phase region above '~1000°C 


PtTc3-D0i9 


Pt-Zr 


two- phase region above 600° C 


PtZr2-C16 (uncertain) 


Pt-Zr 


two-phase region above 600° C 


PtZr3-A15 (uncertain) 


Rh-Ru 


solid solution Ru-A3 


RhRu2'"'°'° 


Rh-Ru 


solid solution Ru-A3 


RhRuP"""*" 


Rli-Tc 


solid solution Tc-A8 


KhTr;!-DOii) 


Rli-Tc 


solid sululioii l c-.V3 


llliTc-ijl9 


Rh-Tc 


two-phase region above 1000° C 


Rh2Tc2-ZrSi2 


Rh-Ti 


two-phase region above 600° C 


Rh2Ti-C37 


Kri- 1 


two- phase region above 0°C 


Kii 1 ( (uncertain ) 


Kn-/jr 


two- phase region above 0°C 


ixnZji4-L) La/ tie-ime 


Rh-Zr 


two- phase region above 0°C 


RhZrs-FCCj^^j^ /tie-line 


Rh-Zr 


two-phase region above 0°C 


Rh2Zr-C37/tie-line 


Ru-Tc 


disorder solution (Ru,Tc)-A3 


RuTc3-D0i9 


Ru-Tc 


disorder solution (Ru,Tc)-A3 


RuTc-B19 


Ru-Tc 


disorder solution (Ru,Tc)-A3 


RU3TC-DO19 


Ru-Ti 


two-phase region above 600° C 


RuTi?J'"°'°, Appendix (XI) 


Ru-Ti 


two-phase region above 600° C 


RuTi2-C49 


Ru-Y 


two-phase region above 0°C 


RUY2-CI6 (uncertain) 


Ru-Zr 


two-phase region above 400° C 


RuZr4-Dla 


Tc-Ti 


disorder /3Ti-A2 


TcTif Appendix (XI) 


Tc-Ti 


(lisordor ;'-?Ti-A2 


TcTi2-C19 


Tc-Ti 


two-phase region 1 cTi^;\^ 


Tc2Ti-(.:il5 


Tc-Y 


no information 


TCY3-DO11 


Tc-Zr 


no information 


TcZr4-Dla 


Tc-Zr 


no information 


TcZr2-C49 



TABLE 7b: Experimental solid solutions, two-phases and "not studied" regions and possible ab initio predictions 
(Table 7a + Table 7b = 96 entries). The table starts in the previous page. 
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D. Experimental compounds that could not be checked by our calculations 



Experimental compound unavailable compound 


o y o (jt,iii 


llllClllidil L KjO LilLi 


TX\J (ill/til/liU LKjOUHj 


A o- A4o- 
rt.g-iVig 




UlililiUWll fcJX|Jt;i ILllKiLHiClil KjiJlLLjJyJUlLH 




Ap-rT'VT A- Ao-n OHi A 
-'T-Sol i 14 -'T-feol vJ*-ll4 


LlllCLVClllClLflC? KJljKJljy 


Au-Cd 




iinkiiowri Gxpcrimditcil compound 


Au-Cd 




un9.v&.ilciblc prototype 


Au-Cd 




11 n k'Ti(~ns7n r>Yr~iori tti oti t a 1 r^i^im ririn nd 

UllJVllWW 11 V L^l llllL^llL CLl L^Wlll LJLJ U.llV_l 


Au-Nb 


A 1 1 n l\r r»'> ( 1 1 n Imrkixrn i 


11 Ti k'Ti(~iw7n r>Yr~iOvi XY^ on t a 1 oi^iTn r~imi nd 

UllJVllWW 11 V L^J. llllL^llL (Xl L^Wlll L/LJ Llllvl 


Au-Ti 


AiiTi-RI Q 


1 in 1 1 n r~»l r> rMTit r»t\/'r~»ri r^fr— ti+rnr'n lomctvi?' R T Q 

Lill<X V (XIlOj UlL^ L/lLJliLJliyL/L^ VJll O LL^l^lllv^lllC J ^ 


Au-Y 


U. X J 1 Li-lliVllLJ W 11 / 


1 inlf"n(~i\s"ri ovr>OT"i rri on t a 1 rrim T^r^n n d 

UlliVllW W 11 t_.jvl_/L^i illlL^llLrdil L'^lllL'V^ LlllLl 


Au-Zr 


Aii^7r ('nP9n - Primal 

-Tl- LL/^ Z-IL X lllllCt f 


1 m atra i 1 a ril T^yoktott'^rT^fi 
uiidiVciiiiciL'ic ^1 \jij\j\jy 


Au-Zr 


AnmZr-? 

-iT-U J^(J iJl y Klj±.tj'-t J 


LI1ICI1VCI1IICILI1C7 ^jx \jij\j\jy 


Au-Zr 


A 11 vi f.vr" / 1 1 n l^T^/^Tsm 'r\T'f^"J"f^"i""\rT>o 1 


LllililiUW 11 tJXpcl lllltJllbdl CUllipwUllU 


Cd-Pd 


'^1 / 1 1 n CX^XTW \ 

y X I uiiiiS-iiL* w 11 J 


U11K.11UW 11 fcJA.pt?l llllt:;lltcll UUllipOUllU 


Cd-Pd 




nnavai 1 a Vtlp ■nrrttotvnp 


Cd-Pt 


r 1 LillJVll*^ W 11 y 


1 1 n 1^'n r^urn f'VTii^T'imfin+al onnrrii^nnd 

LllUvlHJ W 11 LJCl llllC^llUOil Ky\JlLLlJ\J LIllU 


Cd-Pt 


'V'l - 1 1 1 Ti L"Vi ms7"Ti 1 
y^ i LiiiiviiLJ w 11 J 


1 in k'n(~iw7n OYr~iori ni on^"al rrtmnr^nnd 

UllJVllWW 11 VL/L^J. llllL^llUOil L/lJlllL^V^ LlllU 


Cd-Y 


Cd^r Yi 1 -Crl/|r:Smi 1 

v_it_i45 1 11 v_jU4501iiii 


1 in p\\rp\ 1 1 n r»l Timt nt^/'r^o 


Cd-Y 


'^^-'^SS J- 13 ^ LllS^'-'-oS 


1 in scvrpt 1 1 n r\ \ r~imt r»t\/'rio 

UllfX V (XIlOj L/IL, L/lLJliLJliyL/L^ 


Cd-Y 




nn aim i 1 a ril TMTk+rkt^rT^f^ 

LlIl«iVaillClL»lC ^1 KJliKfby \J" 


Mo-Pt 


C l_yL/iy CLl^^JVC XiJLfJ V-i 


nn a'^T'ai 1 a r»lfi T\T'<Tl"iTl"'\7"r\fi r»Tr— cl"i^ir*nir»m<3t7'\7 TlOn <-. 

U.lldVCIillClL'lC ^HJlyVjIj^^C Ull iSlyV^H^lllUlllCUi _y ±^\J\Kf 


Nb-Pt 




1 iTi m^a 1 1 a 111 riTTitntTrio /t T^Tlac;o 


Nb-Hli 


(' l-^' '-'lr>l''lli / f> 


Ill IM\"M1 M 1 ku' 1 ^IV^I / ^l "\"1 


Nb-Rh 




LlllcLVcLlldUlt^ piUlUiypt) 


Nb-Tc 


NhTro Nhn isTm oK-rvMn 

IN IJ _1 ,iN U(J, 15 _L >--(J,S5 CtiVJ-ll 


nn a'^T'ai 1 a r»lfi Tiirkt rkttrTifi 

UlllOi VOillCt L*!" Ml ^ LL* L y 


Pd-Ti 


1 "-i,-} 1- ij ^ '-^z • 


un9.v3yilciblc prototype 


Pd-Ti 


PdrTi -) CI K 


1 in a\ra 1 1 a r^l TiiTtt nt\rr^o 


Pd-Y 


iinlrnrtism at" fifi fivri Pd 

UllxvllLJ W 11 <JjL v^U.V^/U 1 


1 in k'n(~iwm ovriori m on ("a 1 ooim rinn nd 

LillJVllWW 11 V L/L^l llllL^llL cXl L^Wlll L/LJ U.11^ 


Pd-Y 


unknown at ^^7 Pd 

1111J\.1H_J W 11 OjLi OI.^J/U 1 LI 


iink'nortsrn fi"vT^fivimon("al oi^inrnrtnnd 
uiiiviiww 11 L^ci iiiiL^iiLrcxi L,*_jiii u.11^ 


Pd-Y 


Pd 1-1 Y^i- f 1 m Irvi msT^Ti 1 
1 1 5 iLiiiiviiu vv 11 J 


iink'nou^n <3Yn<3TiTvi(3n^"al ooimr^rtiind 

UllJVll^ W 11 CA-LId lllldlUcXl L^Wlll L/LJ Llll^ 


Pd-Y 


PdnVo-ViR 1 ("nnknown^ 


nnk'nouj'n f^vr^tsTi m <^n+ a 1 orinrnokiind 

UllJVllW W 11 \^J\.lJ\^L lllltl^llUOil ^Ulll UlllU. 


Pd-Y 


Pd^YQ-hR14. ('nnknnwnl 


nnk'nouT'n t^vr^fiTimfintal rrinTr^r^nnd 

LilliVll^ W 11 \^J\.lJ\^L lllltZ^llLidil \u\JiXXlJ\J LlllU 


Pd-Y 


/^Pd If n 1 n n Ifn r»TH7"n 1 

l-tX *-*0 !2 I Lllllvll'-' w 11 1 


nnk'nr»TH7'n f^vntsTim^sntal rrimni^nnd 

U.lllvll'-' W 11 C^^JCl lllldlLidil (jUlll^JW UllU 


Pt-Ti 




UlldVclllclUlt:? piOLObypt; Ull oHJlt-lllUlllcbiy ±J±CK 


Pt-Ti 


IT L2 -L 1 rx(_-oi2 


n n a"\7"a 1 1 a ril 1^ f^v^'^'t" r^"t"'\ '■T\t^ 
U,lldVdlldl.Jlt; |.Ji (j bKJ b V JjC 


Pt-Y 


Pt4Y3-Pd4Pu3 


iina,vaiilcible prototype 


Pt-Y 


Pt4Y5-Pu5Rh4 


unavailable prototype 


Pt-Y 


Pt3Y7-D102 


unavailable prototype 


Rh-Ti 


RhsTi (unknown) 


unknown experimental compound 


Rh-Ti 


RhsTia-GcsRhs 


unavailable prototype 


Rh-Y 


RhaY (unknown), hP24-P63 /mmc 


unknown experimental compound 


Rh-Y 


Rh3Y7-D102 


unavailable prototype 


Rh-Y 


RhsYs (unknown) 


unknown experimental compound 


Rh-Y 


Rh2Y3 (unknown), tI140-I4/mcm 


unknown experimental compound 


Rh-Zr 


Rh4Zr3 (unknown) 


unknown experimental compound 


Rh-Zr 


Rh5Zr3-Pd3Pu3 


unavailable prototype 


Ru-Y 


Ru2Y5-C2Mn5 


unavailable prototype 


Ru-Y 


RU25Y44 (unknown) 


unknown experimental compound 


Tc-Ti 


X-A12 


unavailable prototype 


Tc-Zr 


Tc6Zr-A12 


unavailable prototype 



TABLE 8: Experimental compounds that could not be checked by our calculations, because the proper struc- 
ture prototype or concentration is not known or not in our Ubrary. The table contains 21 "unknown experimental 
compounds" and 27 "unavailable prototypes" (48 entries total). 
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E. Experimental compounds in disagreement with other ab initio compounds or two-phase regions 



Experimental compound <^ ab initio wrong compound or two-pliases 


System 


Experimental result 


Ah initio result 


AE, (us-lda) 
(iiic\"/aLom ) 


Ai^Y ('f)aw-(j(ja) 
[nic V /aioni) 


Note 


Au-Nb 


AuNb3-A15 


two-pliase region 


7 


6.5 


(a) 


Au-Y 


AuY-B2 [170] 


AuY-B33 


26 


25 


(b) 


Cd-Nb 


Cd3Nb-Ll2 


immiscible system 


70 


>100 


(c) 


Cd-Pt 


CdPts - a'-Ll2 [181] 


CdPtf Appendix (XI) 


25 


10.4 


(b) 


Nb-Rh 


K(NbRh3)-Ll2 


NbRh3-Al3Pu 


8 


5.3 


(a) 


Nb-Ru 


NbRu3-Ll2 


NbRu3-D024 


8 


2.5 


(a) 


Nb-R,u 


NbRu'-Llo 


two-phase region 


20 


4 


(a) 


Pt-Y 


l'tY-B27 


PtY-B.33 


50 


60 


(c) 


Pl-Zr 


in;iZr.",-I)8« 


t.\vu-])lia.st; it;gitjil 


:',(> 


23 


(c) 



TABLE 9: Experimental compounds (which wc have in the hbrary as prototypes) in disagreement with other ah 
initio compounds or two-phase regions. We inchide the differences between the formation energies of the phases in 
disagreement, both for US-LDA and PAW-GGA potentials. (Gentries). Note (a): the ah initio groimd state is within 
less than lOmeV/atom of the experimental ground state (4). Note (b): the assigned experimental structures are 
poorly justified (2) [170,181]. Note (c): unambiguous significant disagreement (3). 
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XI. APPENDIX: NEW STRUCTURE PROTOTYPES WHICH ARE SUPERSTRUCTURES OF FCC, BCC, HCP 



Systsm 


CdPts, PdPts, 


MoTi, 


MoTis, MosTi, NbsTc, 


RhRu 


RhRu2 




PdsPt 




RuTia, TcTis 






k.^ U kf lOiu ul^^ 


FCC ABq 


BCC AoBo 


BCC ABq 


HCP AtBo 


HCP A2B4 




1 Ivi" ri /If ri /~iTnr\i r* 


\Jl XiLUJL llLllllLliU 


\Jl lillLll llOlllL'lv.' 


' 1 V"! rr^Ti 9 1 
±1 IgOlldl 


\Jl liLUJL LUJLLLUlLi 


Space Group 




imma if 


immm ^ 1 1 


trOLLLL ^i04t 


^_'mcm jfoo 


J^CctloUli oyiuuoi 


oC8 


oI8 


oI8 


hP4 


0C12 


Primitive 












vectors (cart.) 












ai/a 


(1,-1/2,1/2) 


(3/2, 1/2, -1/2) 


(3/2, 1/2,-1/2) 


(1/2, -V3/2,0) 


(-1/2, 3v^/2,0) 


32 /a 


(-1,-1/2,1/2) 


(1/2,3/2,1/2) 


(1/2,3/2,1/2) 


(1/2, v^/2,0) 


(-1/2,-3^3/2,0) 


as /a 


(0-1/2,-1/2) 


(-1/2, -3/2, 1/2) 


(-1/2, -3/2, 1/2) 


(0,0,2^873) 


(0,0,^873) 


Atomic 












positions (tract.) 












Al 


(0,0,0) 


(0,0,0) 


(0,0,0) 


(0,0,0) 


(0,0,0) 


A2 




(1/4,3/4,1/2) 




(1/3,2/3, 1/4) 


(5/9,4/9, 1/2) 


Bl 


(0, 1/2, 1/2) 


(1/2, 1/2,0) 


(1/4,3/4,1/2) 


(0,0,1/2) 


(2/9,7/9,1/2) 


B2 


(1/2,0,1/2) 


(3/4, 1/4, 1/2) 


(1/2,1/2,0) 


(1/3,2/3,3/4) 


(1/3,2/3,0) 


B3 


(1/2, 1/2,0) 




(3/4, 1/4, 1/2) 




(2/3,1/3,0) 


B4 










(8/9, 1/9, 1/2) 



TABLE 10. Geometry of FCC, BCC, HCP superstructures which appear as new prototypes in our study. Positions 
are given as unrelaxed positions in the parent lattice. 

XII. APPENDIX: RELAXED STRUCTURE PROTOTYPES 



System 


MoZrs (AB3) 


MosTi* (AB5) 


MoZrs* (AB5) 


NbsRu* (AB5) 


Lattice 


Orthorhombic 


Monoclinic 


Monoclinic 


Monoclinic 


Space Group 


Imma #74 


C2/m #12 


C2/m #12 


C2/m #12 


Prototype 


MoZrs 


MosTi see note* 


MosTi see note* 


MosTi see note* 


Primitive vectors 

(a, fe,c) (A) 
(a, /3, 7) degrees 


(5.678, 5.678, 5.678) 
(145.5, 128.9, 63.1) 


(5.192,5.192,9.879) 
(139.8, 139.8, 35.0) 


(5.695, 5.695, 10.881) 
(138.5, 138.5,34.1) 


(5.339,5.339, 10.178) 
(139.8,139.8,35.1) 


Atomic positions (tract.) 
Al 
Bl 
B2 
B3 
B4 
B5 


(0,0,0) 
(0.259,0.722,0.537) 
(0.524,0.500,0.024) 
(0.815,0.278,0.537) 


(0,0,0) 
(0.166,0.166,0.833) 
(0.334, 0.334, 0.665) 
(0.500,0.500,0.499) 
(0.667,0.667,0.334) 
(0.834,0.834,0.166) 


(0,0,0) 
(0.205,0.205,0.864) 
(0.349,0.349,0.692) 
(0.495,0.495,0.513) 
(0.642,0.642,0.313) 
(0.800,0.800,0.160) 


(0,0,0) 
(0.170,0.170,0.836) 
(0.335,0.335,0.672) 
(0.501,0.501,0.500) 
(0.666,0.666,0.328) 
(0.831,0.831,0.164) 



TABLE 11. Geometry of relaxed structures which appear as new prototypes in our study. Positions are given as 
fractional positions of the primitive vectors. Note *: MosTi, MoZrs, and NbsRu arc slight distorsion of the prototype 
MosTi'"'"*", monoclinic lattice, space group C2/m #12, with atoms in the following Wyckoff positions: Ti in 2a, 
Mo(l) in 2c, Mo(2) in 4i {x = 1/3, z = 1/6), Mo(3) in 4i {x = 2/3, z = 1/3), and primitive vectors ai/a=(3.171, 1, 0), 
a2/a=(-0.315, 1, 0), a3/a=(-0.638, -0.201, 1.987) [17]. 
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